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Abstract 

The thermodynamical partition function for a gas of color-singlet bags consisting of fundamental 
and adjoint particles in both U(N C ) and SU(N C ) group representations is reviewed in detail. The 
constituent particle species are assumed to satisfy various thermodynamical statistics. The gas 
of bags is probed to study the phase transition for a nuclear matter in the extreme conditions. 
These bags are interpreted as the Hagedorn states and they are the highly excited hadronic states 
which are produced below the phase transition point to the quark-gluon plasma. The hadronic 
density of states has the Gross- Witten critical point and exhibits a third order phase transition 
from a hadronic phase dominated by the discrete low-lying hadronic mass spectrum particles to 
another hadronic phase dominated by the continuous Hagedorn states. The Hagedorn threshold 
production is found just above the highest known experimental discrete low- lying hadronic mass 
spectrum. The subsequent Hagedorn phase undergoes a first order deconfinement phase transition 
to an explosive quark-gluon plasma. The role of the chiral phase transition in the phases of the 
discrete low-lying mass spectrum and the continuous Hagedorn mass spectrum is also considered. 
It is found crucial in the phase transition diagram. Alternate scenarios are briefly discussed for 
the Hagedorn gas undergoes a higher order phase transition through multi-processes of internal 
color-flavor structure modification. 



I. INTRODUCTION 



The thermodynamical description of the strongly interacting hadronic matter can be 
approximated by a free gas with a modified level density of hadronic states. This new level 
density is given by the statistical bootstrap equation. A solution to this equation exists 
only for some range of parameters . It has been argued that the statistical evaluation 
of the mass spectrum in the bag made using the micro-canonical ensemble behaves the 
strong bootstrap model 3|, |J] . Naively, the thermodynamics implied by this hadronic mass- 
spectrum has a limiting temperature T . The bootstrap model with internal symmetry 
of the fireball provides subsidiary variables and allows for new types of phase transitions. 
The fireballs, described by the bags of confined quark and gluon components, stand for 
excited exotic hadrons. They are denoted as the hadronic bubbles or more definitely as 
the Hagedorn states. The Hagedorn states are defined by the fireballs which are consisting 
quarks and gluons in the color-singlet states. The confrontation with the experiment is 
made by projecting the color-singlet states for the non-Abelian color gauge symmetry. In 
the simple model, it is assumed that there are no interactions among the color charges in 
the model except for a color-singlet constraint from the confinement assuming the results 
are stable to the introduction of small interaction terms. Gross and Witten [7J have studied 
the action for fundamental particles in the large- N c limit of the two-dimensional U(N C ) 
lattice gauge theory. Using the functional integral method of Brezin, Itzykson, Parisi and 
Zuber [8] in the large iV c -limit. The action is explicitly evaluated for all fixed weak-strong 
coupling \ ws — N c g 2 by steepest-descent methods. In this limit, particular configuration 
totally dominates the functional integral. The Vandermonde determinant, which appears 
in the measure of the lattice gauge theory also contributes a potential term. Their solution 
suggests the existence of Gross- Witten point for a third order phase transition from weak- 
to strong-coupling \y/s interaction in the large- N c four dimensional lattice gauge theory. 
They have also mentioned that the occurrence of such a phase transition would not mean 
that the large- N c theory does not confine but more precisely that the occurrence of such 
a phase transition are not described by the same analytic functions. This means that the 
behavior of the mass spectral density must be modified at the Gross- Witten point. 

The Gross- Witten solution [7] is also obtained for quarks and gluons by solving the 
singular integral equation. The third order phase transition point is determined when the 
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solution of 
one 



;he singular integration is evaluated over an open contour rather than a closed 



The method given by Azakov, Salomonson and Skagerstam 



111 ] leads to 



the singular integration over either closed or open contour. They have considered gluons 
and quarks by solving the resultant singular integral equation and they have obtained a 
third order phase transition point. They have derived the phase transition point considering 
two possibilities. The first one is the integration over a closed contour while the other 
possibility is the integration over an open contour. The phase transition takes place when 
the integration is performed over an open contour rather than the closed one. These two 
possibilities complicate the solution in the realistic case. The order of the phase transition 
is assumed to be the same even for arbitrary constituent mass It is interested to note 



here that the authors of Ref. 



ll| have shown the existence of critical chemical potential [i c 



such as that for T > the physical properties for the low-lying spectrum are unaffected 
by the chemical potential |/i| < /i c . This means that generalizing and extending the results 
for zero chemical potential to the diluted nuclear matter with a small chemical potential is 
satisfactory though more thorough investigation is required in particular when the chemical 
potential exceeds the critical value. The low temperature phase has free energy of 0(N®), 
interpreted as a gas of mesons and glueballs while the high-temperature phase has free 
energy of order O(N^) which is interpreted as exotic color-singlet bags of adjoint gluons and 
of fundamental and anti-fundamental quarks. 

The deconfining phase transition in SU (N c ) gauge theories have been studied at nonzero 

nnn 

temperature using a matrix model of Polyakov loops |13l . 1141 . 1151 ] . This model has been also 
extended for nonzero density. The effective action for loops starts with a potential term. At 
large N c , the action is dominated by the loop potential. It has been demonstrated that the 
Gross- Witten model represents an ultra-critical point where the deconfining transition at 
N c = oo is close but not at the Gross- Witten point. Although masses vanish at the Gross- 
Witten point, the transition is found of first order and it has been suggested that at finite 
N c , the fluctuations can derive the theory much closer to the Gross- Witten point [l3, 14, 15]. 
The thermodynamics of quark iets with an internal color structure has been considered 

nnn 

in the context of a one dimensional quark gas [161 . 1171 . I18I ]. It is considered based on an 
exact Hamiltonian formalism where the quarks are treated as classical particles but their 



interactions through the group theory U(N C ) or SU (N c ) gauge fields are treated exactly 

H,H HQ- 



is 
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The class of theories on a compact spatial manifold with fundamental flavors and adjoint 
gluons are also found undergo a third order deconfmement phase transition at a temper- 
ature T c proportional to the inverse length scale of the compact manifold {21 1. The same 
thermodynamic behavior of the deconfmement transition for large N c at zero and weakly 
coupling is found for a wide class of toy models such as M = 4 super Yang-Mills theory. 



is argued that AdS black hole thermodynamics is related to Hagedorn thermodynamics 22]. 



Sundborg 



t 



221 ] has calculated the partition sums for M = 4 super Yang-Mills on S 3 and 



he has discussed the connections with gravitational physics. It have been speculated that 
on a dual string theory interpretation that the deconfmement transitions are always asso- 
ciated with black hole formation and furthermore that the intermediate temperature phase 
associated with second order deconfmement transitions would be dual to a string theory in 



a background dominated by a strange new type of stable black hole [23|, [24], [25(. Indeed, 
the Gross- Witten point and the Hagedorn states remain very exciting and rich physics. For 
example, in the language of AdS/CFT duality, the Horowitz- Polchinski point for a small 
black hole should correspond to a Gross- Witten transition. This can be related to the phase 
transition for the small black hole to a big black hole 26]. 

The internal color symmetry of the bound state remains to be of a color-singlet state even 
for finite temperature and chemical potential due to the color confinement. The discrete 
low-lying hadronic mass spectrum is generated by the broken chiral symmetry. The inter- 
nal color structure of the hadronic bag of quarks and gluons for hadronic states above the 
known discrete low-lying hadronic mass spectrum remains in a total color-singlet state. In 
the realistic calculations, the internal color structure is imposed in the partition function for 
the hadronic bubbles 27] . When the temperature reaches the critical one, color is expected 



to be liberated and chiral symmetry is restored. QCD predicts a phase transition from the 
hadronic phase to a deconfined quark-gluon plasma. The hadronic phase consists of the 
whole hadronic mass spectrum including resonances of all known particles. Hardcore repul- 
sive forces can be represented by the excluded volume. The effects from strong interactions 
are included by adding a free gas of hadronic bubbles which are bags of constituent quarks 
and gluons with a specific internal color-flavor structure. The hadronic bubbles remain in 
overall color-singlet states, despite the complexity of the internal color-flavor structure. The 
bubble's size grows up but retains its own internal color symmetry with increasing the bary- 
onic density at low temperature. The volume fluctuations are expected to be suppressed 
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whenever the bags start to overlap with each other at large chemical potential. When the 
temperature increases, the bag's surface is deformed and smeared out until the bubbles 
dissociate eventually at the critical temperature. It is expected that the bubble's volume 
fluctuation increases at high temperature in contrast to the cold dense hadronic matter 
where the bubbles tend to squeeze each other. The strength of the volume fluctuations is 
presumed to modify the details of the phase transition diagram [281 ]. 

The (grand-) canonical ensemble and its Laplace transform to the micro-canonical en- 
semble for gluonic bags or glueballs was derived by Kapusta js, 4] without imposing any 
color constraint. The internal symmetry constraint was originally introduced for the sta- 
tistical bootstrap model 0, E]. The bootstrap density of states can be derived from the 
MIT bag model. In the hadronic phase, the highly excited fireballs derived from the boot- 
strap equation are hadronic bubbles of confined quarks and gluons in a net color-singlet 



state 
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34J. Gorenstein et. al. 
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36, 



37 



38 



have studied the gas of 



bags using the isobaric partition function. Auberson et. al. [39(] have studied the phase 
transition to the deconfined quark-gluon plasma by considering the asymptotic Gaussian 
volume fluctuations. The order and the shape of the phase transition have been studied in 



details 



28] . It is shown that the color-singlet constraint imposed on the bag states is not only 



the critical condition for the existence of the phase transition to an explosive quark-gluon 
plasma. Therefore, the appearance of the deconfined phase transition depends essentially 
on the volume fluctuations beside the internal color structure constraints. Recently, it has 
been suggested that the phenomenology of hadronic bubble internal structure decides the 
order of the phase transition for low chemical potential and high temperature 28J, [40] . 

In the present work the partition function is written as a function of the thermal running 
parameter A. This thermal running parameter is related to the Gross- Witten weak-strong 
coupling X ws as follows \/N% = l/\ws- Therefore, the weakly interaction corresponds the 
large thermal running A > A cr j^j ca j while the strong interaction is analogous to the small 
thermal running A < -^ cr itical - 

The outline of the present paper is as follows. In section [Til we review the partition func- 
tion in the Hilbert space for various statistics for the gas of particles with a specific internal 
color symmetry. The partition function is derived for an ideal gas. The sum of states for 
particles confined in the cavity is presented for the convenience in the realistic calculations. 
The extension to a specific geometry and other degree of freedom in the conformal fields 
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such as the super-symmetry Yang-Mills theory is mentioned. The invariance measure which 
is used to project a specific internal color state is presented in Section IIHI The invariance 
measure is given by the Vandermonde determinant for the U(N C ) and SU(N C ) theoretic 
groups. It is shown that the Vandermonde determinant have two different asymptotic ap- 
proximations. The first asymptotic approximation is appropriate when the color charges are 
distributed uniformly over the entire circle. While the second one is more appropriate when 
in the color charges become more dominated in a specific physical range. In Section IIVt we 
review the spectral density of the color eigenvalues method. It has been studied by Gross 
and Witten to discover the critical Gross- Witten point. This method is reviewed for the 
fundamental particles and the phase transition point is given. In Section |Vj we present an 
alternate method to derive the critical point for the phase transition. This method is based 
on the Gaussian-like saddle points method and it is more appropriate in the theoretic nuclear 
physics and realistic calculations. It has be considered by several authors [28l.l29l.l30l.l34l.l41]. 



In Section rVTl the gas of adjoint particles is reviewed in the context of two different methods. 
Since in the realistic situation we have a gas of fundamental as well as adjoint particles, the 
physics of the phase transition for a gas consisting of different particle species and satis- 
fying various statistics is studied in section IVHI At first, we demonstrate the calculation 
for the Maxwell-Boltzmann's statistics for fundamental and adjoint particles. The realistic 
confined quark and gluon bag is considered in detail. The adjoint gluons which are obeying 
the Bose-Einstein statistics and the fundamental quarks and antiquarks which are obeying 
the Fermi-Dirac statistics are considered in computing the partition function and finding 
the phase transition's critical point. The micro-canonical ensemble as a density of states is 
presented in section IVIIII The micro-canonical ensemble is derived from the inverse Laplace 
transform of the grand-canonical one. The critical mass for emerging the fireballs which 
are obeying the Maxwell-Boltzmann statistics is demonstrated at first and then the realistic 
critical mass for emerging the Hagedorn states is given. In section IIXI the density of states 
is derived using the statistical evaluation of the micro-canonical ensemble. The derivation is 
given along the lines of Kapusta Q| but the internal color symmetry is considered explicitly 
in our calculations. The thermodynamics for the gas of bags with the excluded volume 
repulsion is presented in section [X]. We summarize the conditions for the phase transition in 
the context of the isobaric partition function construction. Furthermore, various scenarios 
for the phase transition are given. The role of the chiral symmetry restoration phase tran- 
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sition in the hadronic matter and quark-gluon plasma phase diagram is given in section IXIl 
Finally, our conclusion is given in section IXII1 



II. THE CANONICAL PARTITION FUNCTION 



The canonical ensemble for the gas of fundamental particles in the Hilbert space is given 
by the tensor product of the Fock spaces of particles and antiparticles 



z q -M 



(1) 



ISU(N C ),U(N C ) 

On the other hand, the canonical ensemble for the gas of fundamental and adjoint particles 
is given by the tensor product of fundamental particles, antiparticles and adjoint particles' 
Fock spaces. This tensor product reads 



Zqqg ((3) 



dfi(g) \Tr q U q (g)e-^"] \Tr q U q (g) e -^] \Tr g U g (g)e-^ 



(2) 



'SU(N C ),U(N C ) 

where the notation Tr ?i q )9 is the trace in each Fock space. We adopt the notation Tr in 
order to distinguish the single particle statistics from the trace tr c over the color degrees of 
freedom. In each Fock space there exists basis that diagonalizes both operators as long as 
H q and U q (g) commute. Let \a,j > be the one-particle states of such a basis, where a labels 
the eigenvalues of H q and j labels those R(<?) of U q (g), then the diagonalized eigenstates 
read 



< oc',j'\H q \a,j > = 5 a ^5jj>E a , 
<a',j'\U g (g)\a,j > = S^SjjrKjj^g). 



(3) 



Any configuration of the system is defined by the set of occupation numbers {n a j}- The 
additivity of H q simply means that 



< {n aJ }\Hq\{n aJ } >= ^2n aJ E a . 



a -J 



Then, using the basis |{n Q j} > to evaluate the trace, we readily obtain 



I'd, J 

1 Boson 



TrU q (g)e-^ = £ II(^^)^)e 

K,j} a J 

+ e n(-i)^ +l Ka?)e-^); 

{n a ,j} a J 



a, 3 

Fermion 



(4) 



n e K(^)e-^) n B ( 



Boson 



+ n e (-ir^ +l (r^k^);: 



Fermion ' ^ 

Hence, the partition function in the each Fock space decomposes to Bosonic and Fermionic 
single-particle functions. It is reduced to 

Trstates^K^ = Tr FD U q (g) e -^ + Tr BE U q (g) e -^, 
-> Tr FD U q (g)e-^, 

-> TV BE C/,( 5 )e-^, (6) 

where the subscripts FD and BE indicate Fermi-Dirac and Bose-Einstein statistics, respec- 
tively. The resultant ensemble decomposes either to Fermionic or Bosonic Fock space or 
even the superposition of Fock spaces with different statistics. The Hilbert space of gas, 
which is consisting of particles satisfying Fermi-Dirac and Bose-Einstein statistics, has the 
structure of tensor product of Fock spaces obeying Fermi-Dirac statistics and another one 
obeying Bose-Einstein statistics. The trace over Fermi-Dirac statistics is given by 



Tr FD U q (g)e-^ = J] [l + R(9)e^ E " 



= exp 



+tr c ]Tln [l + R(g)e-P E ° 



(7) 



while the trace over Bose-Einstein statistics is given by 

1 



Tr BE U q (g)e-^ = ]J 



\?[1-R(g)e-P*-Y 



exp 



-tr c J> [l-R(g)e-P E ° 



(8) 



where tr c is the trace over the color degree of freedom. In the realistic physical situation, 
the single-particle partition function satisfies either Fermi-Dirac or Bose-Einstein statistics. 
The quarks and antiquarks satisfy Fermi-Dirac statistics while gluons satisfy Bose-Einstein 
statistics. The sum of states {jn} for Fermi-Dirac statistics is calculated explicitly as follows 



Tr FD^(#) e 



-/3H q _ exp 



exp 



+tr c ]Tln [l + R(g)e^ E 

a 

+tr c EE 1 (-i) m+1 (R(g)e-^y 

' m. v / 



a m=l 
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cxp 



= exp 



+ E 

m=l 

+ E 

m=l 



-1 



,m+l 



m 



m 



tr c R(9 m )J2 e ~ mf3Ea 



im+l 



tr c R(^)^ F (e- m/3 ) 



(9) 



Tr BE U q (g)e-P"« = exp 
= exp 



In Bose-Einstein statistics, the single particle partition function becomes 

a 

a m =l " b 
_ m =l " l \ a / 



= exp 
= exp 



(10) 



+ £ -tr c R(s"> B (e- m ' 3 ) 

. m =l m 

where - 2 p/B( e ~ /3 ) — S Q e - ^" is the sum of energy states with a specific structure and/or 
geometry either for Fermi- or Bose-particles. Furthermore, in the high energy limit (i.e. 
temperature), the Maxwell-Boltzmann statistics becomes an appropriate approximation un- 
der certain conditions. The single particle partition function in the Maxwell-Boltzmann 
statistics is given by 

= nexp(R^)e-^*), 
= exp ^2 m det [exp (R(g)e~ l3Ea ) 
= exp^tr c lnexp (R(#)e _/3Ea ) , 

tr c R(^)^e-^ 



cxp 



(11) 



The theoretical group such as U (N c ) and SU (N c ) has a special importance in the realistic 
physical application such as quark and gluon bubble. 

The U (N c ) or SU (N c ) internal color group symmetry is introduced in the canonical 
ensemble in order to project the state with a specific internal color structure. The color 
structure of the quark and anti-quark is introduced by the fundamental representation. The 
Fock space single-particle partition function with fundamental representation in the Fermi- 
Dirac statistics reads as follows 



Tr FD U q (g)e 



-P H * = exp 



+tr c £ln [l + R {un (g)e-^ 



9 



exp 



+ 



EE ln [ 



1 + e Mi e 



(12) 



i=l a 



where the trace over group U(N C ) (or SU(N C )) fundamental representation is given by 



tr c Rf un ( 



N c 

E< 



ik&i 



(13) 



The partition function for the fundamental particle and antiparticle Hilbert space becomes 



cxp 



+ E E ln [ x + 2 cxa%)e~ pEa + e~ WEa 

i=l a 



(14) 



On the other hand, the gluons are assumed as gauge particles in the theory and they are 
introduced in the context of adjoint representation. The Fock space single-particle partition 
function for adjoint particles satisfying Bose-Einstein statistics reads 



Tr BE U g (g)e 



cxp 



cxp 



-tr c ^ln [l-R ad j(#)e- 



■0E O 



-3* E E ln I 1 - j^-We-P* 1 *] - N c ^ In [l - e"^- 

i^j a a 



.(15) 



The trace over the adjoint representation becomes 

tr c R adj (g k ) = J2 cos k (6 l -9^ + N c , U(N C ) 



cos k 

i+3 



i,) + (AT c -l), 5[/(AT c ), 



(16) 



for the U(N C ) and SU(N C ) groups, respectively. The density of states for single particle 



has been studied extensively using the multiple reflection expansion method 



42 



43 



cvc 



44] 



The sum runs over the number of states which can fill the one particle states. The sum 



over states 
follows 42 



'or constituent particles those are confined in a spherical cavity is calculated as 



43 



44] 



E, 



p 2 + m + V(r), 
£ = V D J v d 3 r J dEp(E), 

= V D V J dEp(E), 



V 



D 



Vd 3 p 
(2tt) 3 ' 



(17) 



The volume V is the spherical cavity volume where the constituent particles are confined 
inside and it corresponds the total occupational space. The one particle degeneracy is given 
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by T>d = N c T>d where N c is the number of colors and T>d is the particle species degeneracy. 
Usually we set T>d = (2j + 1) as the spin multiplicity where j is the spin quantum number. 

Nonetheless, it is interesting to note here that in the case we have a thermodynamical 
system of particles confined in a compactified space with specific geometry and other degrees 
of freedom, the sum of states becomes nontrivial. It has been noted that in the conformal 
field theories, the partition function for scalars and chiral Fermions on S 3 x R are given 



by |22|, 



23] 



and 



(1 - e-Pf ' 



--3 



(18) 



(19) 



;i- e ^) 3 ' 

Hence, the model analysis for the partition function that is given in the nuclear physics for 
the search for a quark-gluon plasma can be generalized to study the conformal fields such as 
pure Yang- Mills theory and J\f = 4 Super-symmetry Yang-Mills theory. The understanding 
of the partition function with various internal structures can shed the light to a new phase of 
matter beyond the quark-gluon plasma. This might be also useful to construct new models 
to explore new physics such as the searching for a dark matter. 



III. INVARIANCE MEASURE 



The invariance measure is essential in order to project the state with a specific internal 
structure. The invariance Haar measure reads 



dfi(g) 



1 ( 1 \ Nc r 
1 / l \ *° /■- 



n 2sin ' 

m>n 

[] 2 sin 

,m>n 




Y[d6 k , U(N C ), 

k=l 

N c / N c \ 

J] [2ir6CE0i))de k , SU(N C )[20) 
j k=i V i=i I 



for U(N C ) and SU(N C ) group representations, respectively. Haar invariance measure for the 
system with homogeneous and uniform distribution over the Fourier color variables, |^| < 7T, 
is written for SU(N C ) as follows 

JVc-l Nc-l 



',21) 
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The Vandermonde determinant, which appears in the measure contributes to the action 
as an additional potential term. The Vandermonde effective potential term is soft when 
the color eigenvalues in the stationary condition distribute uniformly over the entire circle 
< 7r. However, this will not be the case for the extreme condition in particular when the 
color eigenvalues are distributed in a narrow domain \9i\ « ir. Under this condition the 
Vandermonde effective potential becomes virtually singular and the action must be regulated 
in a proper way in order to remove the Vandermonde determinant divergence. The regulated 
Haar measure in the extreme conditions becomes 



J dfi(g) 



1 / 1 



JVJ V27T 



n 

fc=i 



\ 2 

m 

,m>n 



2n5 \Y,0i )de k . (22) 



\j=i / 

This version for the invariance measure accommodates the canonical ensemble in the hot 
medium. The similar results can be verified for the theoretic group U (N c ) representation. 



IV. FUNDAMENTAL PARTICLES: SPECTRAL DENSITY METHOD 



The canonical ensemble for fundamental particles in the U (N c ) or SU(N C ) representation, 
confined in a spherical cavity and obeying Maxwell-Boltzmann statistics, reads 



-f3E(p,r) 



(23) 



where Vd is the single particle degeneracy due to the degree of freedom which is not specified 
explicitly in the calculation such as the spin multiplicity. The trace over the color degree 
of freedom for particles and their antiparticle partners in the fundamental representation is 
given by 



SfunW = ^ tr ( R fun(^) + Rfnn(^))' 



N c 



— cos#j. 



(24) 



c 1=1 



The sum of states for a gas of particles confined in a spherical cavity is calculated as follows 



A = V, 



d 3 p 
(27p 



-f3E(p,r) 



= v t j 



Vd 3 p 



-f3E(p,r) 



(25) 
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The color-singlet state for the canonical partition function for a gas of fundamental particles 
and antiparticles is projected as follows, 

1 ~(iV c 2 -7V c )/2 TT f dOk 



Z(A > = m 2 



x exp 



fc=i 



1 Nc Nc 

2 E E lnsin2 

n=l m=l 



9n 9 m 



V tr ( R fun(s) + Rfun(2) 



(26) 



The multiple integrations in the partition function are evaluated using the steepest descent 
method. The partition function is reduced to 



N c N c 



EE In 

n=l m=l 



6 n — Or 



sin 



1 



+ 2A— £>S0 r . 



c n=l 



(27) 



The saddle points 9 n are determined by the stationary condition 



2A— sin 9a 

N r 



E cot 



9i — 9 r 



(28) 



The above equation is turned out to be an eigenvalue problem for the large iV c limit and the 
saddle points are determined by the eigenvalues. Hereinafter, we introduce A = N^X. The 
stationary condition becomes 



1 



N c 



2A sin 9i = — cot 
N r . ^ 



9i — 9 r 



(29) 



The solution for the stationary condition that is given by Eq.( l29l depends basically on the 
value of the thermal running parameter A. It can be shown numerically that the solution of 



9j — 9 q 



< 



2tt 
N c - 



9i for A < 2 is distributed uniformly over the entire circle range < 7r and 
However, this solution ceases to exist for the full range \9i\ < n when the thermal running 
parameter A becomes relatively large and runs over the range A > ~. The eigenvalues 9i turn 
to be distributed uniformly only over the incomplete range \9i\ < 9 C < it. This means that 
the solution for A > | becomes problematic and the action must be regulated thoroughly 
when 9 C becomes small. When that action is regulated, another analytical solution may 
emerge and the change in the analytical solution characteristic becomes the responsible for 
the existence of the phase transition. 

The numerical solution for the set of stationary conditions which are given by either 
Eq. (1281) or Eq. (l29p versus the thermal running parameter A for various color numbers N c 
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with group symmetry U(N C ) is displayed in FigJTJ The resemblance between 't Hooft cou- 
pling Xw-s = Qym^c and the thermal running parameter A = \/N% is given by the rela- 
tion gyM^c = 1/A. The weak 't Hooft coupling corresponds the high-lying thermal run- 
ning parameter A > A cr ^.j ca j while the strong 't Hooft coupling is given by the low-lying 
A < ^critical' ^ ne Vandermonde determinant contributes to the effective action. Con- 
sequently, the saddle point locations are computed using the stationary conditions where 
the Vandermonde determinant is included explicitly in the action as an additional potential 
term. It is shown numerically that the saddle points are distributed uniformly over the inter- 
val — 7r < $i < it for low-lying thermal running parameter A (i.e. strong 't Hooft coupling). 
The characteristic distribution of the saddle points is found uniform for the small A even for 
a small number of colors N c = 2. This distribution becomes more dominant over a narrower 
range \6i\ < 6q < 7r as A increases and approaches the critical value. However, the numerical 
calculation shows that the saddle points cease to exist when the thermal running parameter 
A exceeds a critical value just above A > \ (i.e. weak 't Hooft coupling). For example, in 
the case of N c = 3, the saddle points cease to exist when the thermal running parameter 
exceeds the value A > -^p. Furthermore, it is found that when the number of colors N c 
increases the saddle points cease to exist at smaller value of A. For example, in the case 
N c = 5, the solution for these points does not exist beyond the critical point A > ^ip. In 
the limit of large N c , the solution for saddle points breaks down just above the point A = |. 
The Vandermonde term diverges when the color saddle points become more dominant in a 
tiny domain around the origin. This divergence breaks the solution badly when these saddle 
points congregate and converge to the origin. It is found that the saddle points cease to 
exist when A > A cr ^.j ca j. This implies the stationary conditions must break down and the 
action must be regulated in a proper way and a new set of conditions must be re-formed. 
This proves the existence of the phase transition in the system. 

In the large- N c limit, Gross and Witten have solved these equations in the case of two- 
dimensional U(N C ) lattice gauge theory following Brezin et. el. 
tral density of eigenvalues. It is shown by Brezin, Itzykson and Parisi and Zuber 
functional integrals in the large- N c limit can be calculated by steepest-descent methods. In 
this limit, a particular configuration totally dominates the functional integral. The station- 
ary equations in this limit can be replaced by their continuum version by introducing a no 



by introducing^ the spec- 
that the 
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decreasing function 

0n = 0(X) 



n 

n = l,---,N c . (30) 
Hence, the stationary condition becomes 

2Asin^ ) =PjT 1 dyoot ^ ( " 0) 2 " %) ) , (31) 
where P refers to the principle part of the integral. By introducing the density of states 

m = % > 0, (32) 

and the constraint 



dx = / d9p{9) = 1, (33) 

J-0c 

where \9 C \ < it. Hence the stationary condition becomes 

2Asin# = P / d0'p{6') cot — — . (34) 



A. The highly thermal excited matter: > \ 

The solution of density of eigenvalues in the large N c limit in U(N C ) representation has 
attracted much attention, recently. Although the solution for the spectral density of eigen- 
values is difficult in the realistic situation, the approximate solution for this density has 
been derived for some simple specific cases for fundamental and adjoint representations in 
the limit of large N c . In the following, we shall review the partition function using the given 
solution of the spectral density of color eigenvalues. Then in section |V] we introduce another 
approximation to evaluate the same partition function. The other method is to approximate 
the resultant integral to the Gaussian-like integral and this approximation sound more effi- 
cient when more complicate situations are involved. The two methods have been shown to 
reproduce the same results for the large thermal running parameter A > Ao and A = X/N%. 

The first method to evaluate the canonical ensemble is based essentially on the stationary 
equation solution for the density of color eigenvalues and it is written for A > | as follows 



2 ~ 9 
p{9) = -Acos- 

7T 2 



• 2 ( Qc\ ■ 2[® 

sin — — sin - 
\2 \2 



n l/2 

(35) 
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where 

sin2 (I) = k £ 1 ( X/N < = x ^) ■ (36) 

The color eigenvalues are not distributed uniformly over the entire color circle range — tt to ir 
but instead, in this solution they are distributed only over a narrow interval \6{\ < 9 C where 
\9 C \ < it. The canonical partition function after weighting the color density of eigenvalues 
becomes 

where the pre-exponential constant is given by C = 2( - Nc ^ N ^ )/2 m Under this stationary ap- 
proximation, the canonical ensemble function is calculated as follows 

,~ oN AT? „„ 3 



Z(X) = exp 



2NfX - ^ In (2AiV c 2 ) + ^ In iV c 2 - °-N> 



■Vi 2 e( 2 ^-^), where A > -. (38) 
2\) ~ 2 K J 



Hence the asymptotic large running thermal coupling solution reads, 

'NT) 2 e (2A-§*2) 



Z{\) - ^ eV»-W. (39) 



B. The diluted and relatively cold matter: < ~ 

At the low temperature, the stationary condition for A < | produces the following density 
of color eigenvalues 

p(Q) = — fl + 2A cos 9), (-7i<6<7i). (40) 
2n ^ ' 

The chemical potential for the total particle number due to the rotated U(1)b symmetry is 
not considered here and it is left for the forthcoming work. In the large N c limit, the system 
is assumed highly compressed due to the large number of colors but this is not the case for 
a finite number of colors such as QCD. In this case, the system is identified as a compressed 
one when the chemical potential for particles numbers becomes relatively large. The density 
of eigenvalues given by Eq. (jlDj) for < | produces the following partition function 

Z{\) = e N ^\ 

= e x2 ' N *. (41) 
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The occurrence of such a phase transition from the diluted and relatively cold matter (i.e. 
low-lying phase) to the highly thermal excited matter (i.e. high-lying phase) would not 
mean that the large- N c theory does not confine. The highly excited thermal matter can be 
interpreted as an exotic hadronic phase dominated by the Hagedorn states. In the lattice 
theory, this would imply that the weak- and strong-Av^s coupling is not described by the 
same analytic functions [7]]. The weak- and strong-Ajy-s are analogous to large and small 
thermal running parameter A/iV|, respectively. For the finite N c , one could continue from the 
strong to the weak coupling. However, one would expect to see a sign of phase transition for 
large enough N c (i.e. in the limit N c = oo) whose manifestation would be a sharp transition 
at \w-s ~ -^w-s(critical) fr° m the weak-coupling to the strong-coupling behavior. As in 
the lattice theory, in the thermal and dense QCD, the strong- and weak-coupling transition 
is analogous to the phase transition from the discrete low-lying mass spectrum particles to 
the highly excited and massive Hagedorn states (i.e. continuous high- lying mass spectrum). 
The Hagedorn states are the color-singlet (confined) hadronic states. This phase should not 
be interpreted as an immediate deconfined phase. The critical Gross- Witten point is the 
threshold point for the Hagedorn states to emerge in the system. The deconfined phase 
transition can either take place immediately when unstable Hagedorn states are produced 
in the system or as a subsequent process when the metastable Hagedorn phase undergoes a 
true deconfined phase transition. 

V. FUNDAMENTAL PARTICLES: GAUSSIAN-LIKE SADDLE POINTS 
METHOD 

We introduce an alternate novel method to calculate the partition function and to find 
the critical point for the phase transition. This method accommodates the sophisticated 
physical problems such as the internal color structure for a gas consisting particles with 
various statistics species in a specific space boundary and with more complicate color-flavor 
correlations and chiral symmetries. It may be useful to find the equation of state for physics 
beyond QCD. 
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A. The highly thermal excited matter: > \ 



In the realistic physical situation, more calculations are involved. It becomes more diffi- 
cult to calculate the density of color eigenvalues under certain conditions in particular when 
the color eigenvalues populate a tiny interval |6> c | <C 1 under extreme conditions. The canon- 
ical ensemble has a special importance for the thermal running parameter A over the range 
A > \ and it corresponds the hot hadronic matter. It is relevant to the relativistic heavy ion 
collisions in RHIC and LHC, and the cosmology such as the physics of the early universe. 

In following, we shall demonstrate a method how to calculate the partition function for 
the coupling parameter ^ > A (in our case A = -^critical = integral over the color 

invariance measure is approximated to the Gaussian-like multi-integrals around the saddle 
points. The color saddle points are dominant in the narrow interval \Q\ < 9 C around the 
origin where \9 C \ <C 1. As far the saddle points are populated in a small domain we can safely 
approximate all the saddle points to be located near the origin. The Vandermonde potential 
does not appear in the action but instead the invariance measure appears as a pre-factor 
power function of the exponential in the Gaussian integral. The integration of the resultant 
multi-Gaussian integrals is straightforward. We shall call this method as the Gaussian-like 
saddle points method. In this limit the invariance Haar measure is approximated to, 



V c) Saddle - II J Mkjfi 



and 



Ln c) ^saddle 



k=i ' 



n>m 



n (o n 

n>m 



Or, 



ISU(N C ) ... k=r 

for U (N c ) and SU (N c ) symmetries, respectively. The partition ensemble is given by 

A 



(42) 



(43) 



Z(X) = J dfji(g) exp 



Vc tr c(RfunG?)+ R fun 

f ( X Nc \ 

J dfi(g) exp 1^2— J2 cos °nj ■ 



(44) 



When the thermal running parameter becomes a relatively large and runs over the range 
A > \//)minimum> ^ ne partition ensemble is approximated around the saddle points which 
become dominant in a small range as follows 

Z(X) » Z (II) (X) (@ A runs over the range A > -^(//)minimum) ' 
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N c 



z (ii)W = / d» (9) saddle ex P 2NcX cos ^ > 



2 TV? A 



^(^)saddle ex P 



iVc 



n=l 

. z n=l . 



(45) 



where A = N%\. In the case of U(N C ) symmetry, the partition function becomes 



%/)(A) = ^F7 ( ^ 



! / 1 ^JV C e 2Af c 2 A 



(2iV c A)' 



n 



T \k=l 



■I'D 



ft (x n -zJ 2 V§£»=i< (46) 



Vn>m 



The Gaussian-like multi-integrals are evaluated using the following standard formula 



< N c 

n 

\/c=i 



\n>m 



exp 



Z n=l 



(27r)^ /2 n nl 



(47) 



n=l 



Hence, the partition function is evaluated as follows 

iV 2 -l 

%/)(A) = 



,2AT?A 



. 2AU . 



1 1 1 



N c 



A > A 



(i/jminimum 



(48) 



and 



Z(ii)(X) 



n 



n l 



A > A 



(//)miiiimum> 



(49) 



2 iy c : n=l 

for SU(N C ) and U(N C ) symmetries, respectively. In order to compare the results of the 
present method with those derived using the density of eigenvalues method (i.e. the spectral 
density method) given by Eq. (l39l . the partition function for U(N C ) is written as follows 



z {II) (\) = c 



2A 



„2Ni\ 



A > A 



(//)mimmum' 



where the pre-factor constant is given by 



N c 

n 



7?! 



(50) 



(51) 



Using the Stirling's approximation Inn! ~ nlnn — 77 + | In 77, + | ln27r, the pre-factor given 
by Eq.l jSTl) is simplified to 



\nC Nc = — ^lniV c 



2 



N C (N C + I) 



- In 77 + n In 77, 

2 n= i n =i 



\nN r . 



N C (N C + 1) 1 rNc 



N c 



dn In 77 + / dnn In n, 



-Nt--N c \nN c . 



(52) 
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The above expression simplifies the partition function to 

@ A > -\//)minimum- (53) 
The above solution is the approximate solution over the range 



/ i 

hm %r)(A) « -~ exp 
Ar c ^large V2A/ 



^2iV c 2 A - ^iV c 2 
t 4 



A - -^critical - A (//)minimum' where -^critical - ^' ( 54 ) 

It is in agreement with Eq. (l39l) . This means that the Gaussian-like saddle points method 
is consistent with the density of eigenvalues method developed by Brezin et. el. [8] for the 
energy range jfe > §. 



B. The diluted and relatively cold matter: < ^ 



The Fourier color angles are distributed uniformly over the entire color circle range — 7r < 

A. < I 

N 2 - 2 



< 7i in the low-lying energy domain -A < as follows 



Z(A) Ri Z^){\) (@ A runs over the range A < A ( - / - )max j mum ) , 
Z { i)(\) = jT d/iQ?) exp [A^tr (R(p) + R* (<?))] . (55) 

The partition function is expanded with respect to orthogonal bases over < ^(/)maximum 
and it reads 

Z {1) {\) = f dfi(g) exp [iV c Atr (R(g) + R*(g))} , 

J ~ 7T L J 

/ w OO OO 1 .1 

M9) E E ^ (NcXtrRig)) 1 - (N c \trR\g)) J . (56) 

w 

/TV / 1 \ ^ O ' 

_ M9) E ( -j) {NS) 1 (trR(^)) J (trRt(^ 

OO 1 

1=1 L - 

e (^ 2 A 2 ) ; 

(over the range A < A critical < A (/)maximum ) , (57) 



i=ij=i"- J- 
Using the orthogonal relations over the U(N C ) bases, the partition function is reduced to 

1 \ 2 / „ r\ 2i 

oo 

1=1 4 - 



where A critica j - | and A (/)max j mum < I. 
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The canonical ensemble with Maxwell-Boltzmann statistics single-particle partition func- 
tion reads 

\ 2" 



cxp 



cxp 



exp 



c 

(2j + l)V s V 

7T 4 iV c 2 /3 6 

(2j + 1)V 5 2 3 T 6 

7T 4 iV 2 



V s3 d 3 p 



-0E 



(58) 



for an ideal gas of fundamental particles confined in a spherical cavity in the three dimen- 
sional space where V$3 the cavity volume in the three dimensional space. The free energy 
approaches the same value as in the free theory for low temperature, whereas we have a 
mesonic-like and baryonic-like matter in the dilute hadronic phase. The canonical ensemble 
in this range does not produce an explosive phase transition. The physics for the two phases 
will be discussed in rather details when we present the micro-canonical ensemble as a density 
of states in section IVIIIl 

We display the function lnZ(A) versus for the gas of fundamental particles in the 
color-singlet state with various numbers of colors N c in Figs. ([2]) and ([3]). The logarithm of 
the canonical ensemble in the domain > | is displayed in Fig. (J2D- The highly thermal 
excited partition parameter ^ > ^ corresponds the weak 't Hooft coupling and it describes 
a hot gas. The partition function in the large A regime is evaluated exactly by the exact 
numerical integration over the Fourier colors over the entire circle range (— 7r < 9{ < 7r). The 
results for the numerical integration are compared with the solution for the spectral density 
of eigenvalues method which is introduced by Brezin et. el. s| as well as the solution of 
the Gaussian-like saddle points method. The two methods are found producing the same 
results and they fit the exact numerical calculations when the thermal partition parameter 



exceeds the critical value > \. In the regime below the critical point Ac 



the two 



methods deviate from the exact numerical one and this deviation increases as decreases 
down below the critical thermal running parameter j^.. The largest discrepancy is found for 
the small number of colors N c = 2. Despite the spectral density of color eigenvalues method 
is originally derived for the large N c — > oo limit, the solution is found satisfactory for a finite 
number of colors even for a relatively small one N c = 2. Furthermore, the Gaussian-like 
saddle points method is found to produce precisely the asymptotic solution for the large 
iV c — » oo limit although it is originally derived for a finite number of colors. However, the 
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extrapolations of both methods to < jfe evidently fail to reproduce the exact numerical 
results in the low-lying energy regime in particular when becomes sufficiently small 
■j§2 "C |. The failure of the large energy extrapolation to reproduce the correct asymptotic 
low-lying energy solution indicates the existence of a possible phase transition near the 
point ^2 — |- This deviation from the exact numerical solution increases significantly 
as decreases below the critical one ffi. The low-lying energy solution for the small 
thermal running parameter domain < | with various color numbers N c =2, 3, 4 and 
5 is displayed in Fig.(j3J). The exact numerical solution is compared with the asymptotic 
analytical solution hxZ — \ 2 /N% for the low- lying energy solution which is derived for the 
small thermal running parameter -^j. The asymptotic analytical low- lying energy solution 
matches the exact numerical one in the range < |. This agreement between the low- 
lying energy solution and the exact numerical one is found satisfactory even for a small color 
number N c = 2. The critical point A cr j^j ca j = -^critical/^c f° r the phase transition from the 
asymptotic low-lying energy solution for the range < ^ to the asymptotic high-lying 
energy solution for the range > is determined when the both asymptotic solutions 
match each other from left and right sides, respectively. Roughly speaking, they split only 
with a small redundant constant due to several approximations considered in the derivation. 
This additive constant vanishes smoothly as the rough approximations are eliminated and 
the solution approaches the exact one. The solutions for the asymptotic small- and large--^ 
energy domains are displayed in Fig. (J4]) . The both solutions are displayed with various color 
numbers N c and they are compared with the exact numerical one. The phase transition is 
found a third order one. This means that the asymptotic low-lying and large-lying energy 
solutions and their first and second derivatives are equal. This implies that both of the 
solutions smoothly tangent to each other and lie above each other. They look like that they 
do not really tend to be intersected but instead they look as they are approaching each 
other. The critical point is the mid-point in the interpolation between the low-lying and 
high-lying energy solutions. The matching of the asymptotic large solution that is 
calculated using the Gaussian-like saddle points method with the asymptotic small A solution 
is displayed in Fig. ([Ha). It is found that at the critical point, the low-lying and high-lying 
energy solutions are split to the minimum with a small redundant additive constant. This 
redundant constant emerges due to consequent approximations adopted in the analytical 
solutions with a finite N c . Although the two solutions deviate from each other by a small 
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redundant constant, evidently, it is always possible to find a smooth interpolation between 
the two solutions. Since the order of phase transition is a third order then the tangent 
slope and its derivative are interpolated smoothly The exact numerical solution is shown to 
match the low- lying energy solution and then when exceeds the critical point, the exact 
numerical solution deviates the asymptotic low-lying energy solution. The exact numerical 
solution continues to match the asymptotic high-lying energy solution instead of the low- 
lying energy solution. This tiny redundant constant should not worry us at all and the 
location of the critical point is realized in the midway of the interpolation between the 
two asymptotic solutions and this constant should be vanished in order to preserve the 
continuation. In this case, the interpolation between the two solutions is smooth and has a 
soft deflection characteristic at the critical point. On the other hand, the high-lying energy 
solution that is obtained by the spectral density of eigenvalues solution for the range > jjj? 
is displayed in Fig.(jHb). It is shown that the low- and high-lying energy solutions smoothly 
tangent to each other and located adjacent to each other for a finite range where the Gross- 
Witten point = \ is located in the midway. The tangent slope and its derivative for the 
both left- and right- solutions are equal. Hence, the both low- lying and high-lying energy 
solutions are scaled to be equal at the Gross- Witten point. 



C. The phase transition critical point 

The critical point for the phase transition is the solution of the spectral density of eigen- 
values method and it comes naturally when the integration over the Fourier color variable 
does not complete the entire circle range but rather covers an incomplete circle range \Qi\ <9 C 
where \8 C \ < 7T. The spectral density method produces a solution for the small thermal run- 
ning parameter -^j < \ and another solution for the large parameter > |. The critical 
point is determined by matching the two solutions. It is interesting to point out that the 
small- and large-A parameters correspond the strong- and weak-couplings in Gross- Witten 
scenario 0], respectively. 

On the other hand, the Gaussian-like saddle points method produces the high-lying en- 
ergy solution only for the thermal running parameter that runs over the interval > |. 
Hereinafter, this is labeled as the solution (II). In solution (II), the Fourier color variables 9i 
are assumed to be dominated in a narrow domain. However, the low-lying energy solution 
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for < ^ is known to be distributed uniformly over the entire circle with respect to the 
color angle. Under this condition, the partition function is evaluated trivially by expanding 
the integral with respect to the group basis powers and using the orthogonality and then 
finally re-summing the resultant terms. Hereinafter, we label the low-lying energy solution 
as the solution (I). The extrapolation of both solutions fails to reproduce the asymptotic 
results far away outside their domains. The point of the phase transition is determined by 
examining the continuity when both solutions match each other along the axis of the thermal 
running parameter A = The thorough study of the solution (II) sheds more information 
about the location of the critical point for the phase transition. The solution (II) is concave 
up for the small parameter The extremum left-hand point for the solution (II), namely, 



i™ J nZ Ni {X) ) = -^n2A + 2A-^,A> A (// )minimum (59) 



is determined by 



d\ Uc-oc W 



— — > A( 7/ ) ex t remum — . (60) 



A (//)extremum 



This point is the bottom of the concave up solution (II): 



d 2 ( Um In^WA)' 



«9A 2 U<=-oo jV 2 



— > + at A^^ m j n j mum — . (61) 

A (i/)extremum 



The left-hand extremum point is found the minimum of the solution and the physical solution 
Z{n){\) is found for A runs over the range oo > A > A^ m j n j mum . This means that the 
solution (II) is the satisfactory solution over the range A > A ( - // - )m j n j mum . The critical point 
■^critical * s l oca ^ e d somewhere above A ( - // - )m j n j mum , 

A > A cr j^j ca j > A ( - // ^ m j n j mum , where A ( . // ^ m j n j mum = -. (62) 

Therefore, when the low-lying energy solution Z(X) w Z^(\) over the range A < -^critical ^ s 
not known, the Gaussian-like saddle points method still has the capability to shed much in- 
formation about the critical point location and the continuous high-lying energy production 
threshold. 

In order to determine the location of the critical point precisely, we need the solution (I) 

for the energy domain A < Am . = 1. The critical point is determined whenever 

v 'maximum 
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the solution (II) from above and the solution (I) from below match each other somewhere 
in the interval that is bounded by the interval 



\//)maximum — ^critical — \/)minimunr (63) 
In our case, the solution (I) reads 

^ " = ^ 2 for ^ < V)maximum- ( 64 ) 

Hence, the critical point for the phase transition between the two solutions is computed by 
matching the extrapolation of the left- and right-hand solutions as follows 

mZ (J) (A) _ \nZ {II) (\) ~_~ 

jp ~ Jp Viticah V° b > 



It is found that -^critical = \~ 



D. Saddle points approximation with the Vandermonde potential 



In order to understand the validity of Gaussian-like approximation, it is essential to show 
that the color saddle points are distributed in a narrow domain (\9 C \ <C ir) as the parameter 
A runs over the range A > -^ cr itical" ^ ne canon i ca l partition function for the fundamental 
particles obeying U (N c ) can be written as follows, 



z(s) = ^n/^«p 



N c N c 



JEE lnsin2 

n=l m=l 



+ 2N C X J2 COS #n 



n=l 



(66) 



The above partition function is extremely difficult to be evaluated exactly. Fortunately, 
it can be integrated over the color variables using the saddle points approximation. The 
exponent under the integral can be written as follows, 



\nZ(\,{6}) 



N c N c 



- ^ ^ In sin^ 



n=l m=l 



n=l 



+ constant. (67) 



In order to avoid the singularity in the exponent, the color saddle points are distributed 
uniformly over the complete circle and do not approach each other (9i ^ 6 m ). However, 
when these points approach each other in the limit 6>j — > 9 m the exponent blows up and 
diverges and subsequently the saddle points integral approximation is badly broken down 
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and another analytic solution must emerge beyond this critical point. The color saddle 
points are evaluated using the stationary conditions 







w .^z(x,{9}) 



o. 



1,2, ...jV c (oriV c -l) 



(68) 



for U(N C ) or SU(N C ). Hence the saddle points are determined by the set of equations 

_ N c (Q _Q \ _ _ 

2N C X sin 9i = £ cot ^ , 9 X ± 9 m . 

n^i V Z J 

These points are satisfying the set of non-linear equations as follows 



(69) 



Mjft^i] = ^ EE cot 

V v c i= i ) iv c i n 



(70) 



This non-linear set is the root of the spectral density that is derived using the spectral 
density of eigenvalues method in the large N c limit. The Gaussian-like integral around the 
saddle points is approximated as follows 

d8i 



Z A 



n 



,r 2vr 



exp lnZ(A,{#}) 



z(\{e})U 



2tt' 



where the exponential elements are determined at the saddle points 



1,2, 



follows 



2 



\nZ(X,{9}) 



Eq. (!7T!) is evaluated around the saddle points using the Gaussian integral as follows 



(71) 

N c ) as 

(72) 
(73) 



Z(\) = A/" exp [lnZ(A,{#})J , 

where the pre-exponential normalization satisfies the integral / ^ = 1 and it is determined 
by 



M = (2tt) 



-N c /2 



det A{9i, 8j 



-1/2 



< 1. 



(74) 



However, it is not always possible to find real values for the saddle points distributed uni- 
formly over the entire circle range for the parameter A that runs along the real axis. These 
real saddle points cease to exist when the thermal running parameter A reaches the critical 
value. The Vandermonde potential characteristics is modified at the critical point. This 
point is the threshold for another solution with a different behavior. 
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VI. ADJOINT PARTICLES 



The color-singlet bags consisting of a gas of fundamental and anti-fundamental particles 
represent the mesonic and baryonic states with no gluonic component. On the other hand, 
the color-singlet bags consisting only of the adjoint particles represent glueballs states with 
no quark component. The low-lying hadronic states are likely mesons, baryons and glueballs. 
The Hagedorn states with fundamental and adjoint constituent particles are essential to 
understand the highly thermal excited hadronic states near and just below the deconfinement 
phase transition diagram. 

The canonical ensemble for the bags consisting only of the adjoint particles projected in 
the color-singlet state reads 

d 3 p 



Z{(3) = J dfi(g)exp 



d 6 r 



v 



(2tt) 



^adjW 6 



-0E(p,r) 



(75) 



where V g is the degeneracy such as the spin multiplicity V g = (2j + 1). The internal color 
structure for the adjoint constituents in the Fock space is given by 

1 



^adjW 



dim, 



-tr 



[ R adj(^)] > 



(76) 



where dim g = N% and (Af? — 1) for U(N C ) and SU(N C ), respectively. The partition function 



for the gas of free particles occupying the volume V is approximated to 

d 3 p 



Xr, 



V g j v d*r 



(2tt)< 



-j3E{p,r) 



(77) 



where E(p,r) is the energy for each constituent particle species. In order to simplify our 
notation, we define the following thermal running parameter 



dim„ 



Hence, the canonical ensemble is simplified as follows 



Z(X g ) = J dn(g)exp 
The adjoint group representation reads 



\ g tr c 



R adjO?) 



(78) 
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R adj (/)] = tr [R fund (^)] tr [R| und (/) 



-1, 



N c N c 

= (N c -1) + Y: E ^sk(9 n -9 m )\ n ^ m , 

n=l m=l 



f N c -1 \ 

V i=l / 



(80) 
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and 



tr 



R adj(^) 



R fund(^ fc ) 



(81) 



: tr [ R fund(^)J tr 

for the groups SU(N C ) and U(N C ), respectively. 

We use the Gaussian-like saddle points method developed in the previous section in order 
to calculate the asymptotic high-lying energy solution for the canonical ensemble with large 
thermal running parameter X g > r j+j i- I 11 this range, the Fourier color variables are 
assumed to be dominant only in a narrow range. The adjoint action near the saddle points 
is approximated to 
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The above approximation reduces the canonical ensemble to 
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for SU (N c ) group representation where the following relation has been adopted 
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The same procedure can be followed for U (N c ) symmetry. The canonical ensemble is reduced 
to 
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Using the Stirling's approximation as done in Eq. (1521) for the large N c limit, the partition 
function becomes, 



W \2\ g , 



N 



2A, 



(@ A 9 runs over the range X g > X g 

(/7")minimum) > (^) 

for the case of U(N C ) representation. The additional constant C ac lj insures that the term 
In Z(u}(\g) /N^ is non-negative quantity and it should satisfy the boundary near the critical 
point and furthermore it throws away any redundant constant. This additional constant 
stems from the normalization transformation in the Gaussian-like integration f°° ^-e° — > 
J-k = 1- Hence, the solution for the large parameter \ g > A 9 ^ // - )m j n j mum reads 



lim lnZ ( n )( X 9 ) 



,jv c 2 ^large 



-iln2A 9 + A 5 -^ + C adj , 



U(N C ) 

> 0, (87) 



in the large iV c -limit. This function is concave up. The value A s( - // - )m j n j mum is calculated 
as follows 

W-4arge<9A 9 JV? L ° ^™ um 2' 

9 (//)minimum 

\iv|-large <m| 7V C / - 

9 (/7)mimmum 

This means that our solution is valid only when \ g runs over the range \ g > A g ( // ^ m j n j mum . 
The value C ac ij is determined as follows 

( , im ^W)>0, (89) 
\7V|^large N£ J 

for the point A 9 ^ // ^ m j n j mum in order to ensure that our solution is non negative in the entire 
A range and also in order to preserve the analytic continuation of the solution. This leads 
to C ac ij = \- Hence, the canonical partition in the range \ g > A 9 (//) m i n i mum in the limit 
large N c limit reads 



1 'NfX g -4 
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It will be shown below (e.g. in Eq. fllOOjl ) that for X g < ^ g ^ ma ximnm wnere 



\(/)maximum < 1> ^ ne canonical ensemble becomes 



lim lnZ(A s ) = lim hxZm(X g ), 
iV|-> large ivj-Jarge 
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(@ X g runs over the rangeA 9 < A 9(/)maximum ) • (91) 
We are in the position to compute the Gross- Witten critical point 

\(/)maximum — \critical — \(//)minimum- (92) 

The critical point for the phase transition is determined as follows 

taZ^A) InZ^X, (o po . nt - = - } 

Ar c 2 ^large JV C 7v c 2 ^large iV c 

It is found that the critical point is located at A gCr ^j ca j = |. This procedure is very useful 
to determine the point for the phase transition for a finite number of colors N c and it 
will be more appropriate in the case of more complicated physical situations are involved. 
Furthermore, it sounds that the Gaussian-like saddle points method is a straightforward and 
easier than the spectral density method. 

The solution of the canonical ensemble for the adjoint particles in the range X g < 
\(/)maximum * s nontrivial because of the color structure for SU (N c ) and U (N c ) group rep- 
resentations. Fortunately, the adjoint color structure can be simplified to fundamental- like 
structure by introducing the Lagrange multiplier. The inclusion of the Lagrange multiplier 
trick reduces the partition function with U (N c ) group structure to 

Z{X g ) = Jdfi(g)exp ( A 9^^ trR adj(^)) > 

= J ^(^)exp(A g trR fun (^)trRJ un (^)) , 
]\f 2 poo iVcS 2 f \ N E / 

= 2X I ^ ^ J M9) 6XP [ 2 ( trR fun(3) + trR fun(^) 

The above equation is decomposed into two parts 

Z(X g )=Z 1 (X g ) + Z 2 (X g ), (95) 
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where 
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(96) 



and 



Z 2 (X 9 ) = §- / d^e-^(0~ N ' /2 e^-^\ (97) 
2X g J£o 

The integration over the color-singlet state is taken into account the splitting into small- 
and large-domains < £ < £o and £ < £ < °°; respectively, where £o = 1- The integration 
over the range (0 < £ < £o), reduces Eq.(l96l) to 
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where A 9 is assumed to run only over the range X g < A g ( J ^ max j mum < 1. The extrapolation 
of the integration upper limit £ — > oo approximates the partition function to 
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This integration is trivial and is reduced to 
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The integration over £ < £ < oo where £ = 1 for the case X g < | < A ff ^ max j mum 
approximated to 
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This integral complements the integral given by Eq. f!99|) which runs over the interval < 
£ < £o- In Eq.( |99|) when £ is extended and extrapolated to oo in term Zi(X g ), then the 
second term is suppressed. 
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On the other hand, the analytical solution behaves differently in the energy domain 
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The integration over the range £o < £ < 00 an d A„ > 1 approximates Eq. fl97|) to 



Z(\ g ) « Z 2 (X 



9)1 

N 2 roc 
c 



2Xg MO 



So 



/«) 



(103) 



where 
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When the saddle point becomes £ sa ddle — the a bove equation is approximated to 
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The saddle point can be found by the following solution 
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The saddle point solution reads, 
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The saddle point can, alternatively, be derived from Eq. (jl06j) as follows, 
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This result is analogous to Eq. (11 181) . The above approximation softens and extrapolates the 
constraint to 



f saddle = 2 V ( with constraint \ g > -, and Saddle ^ 1 

Hence, by substituting the saddle point into Eq. (11041) . we get 
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where we have introduced Sq = V^saddle ^ or our conven i ence - 111 terms of Sq, the partition 
function is reduced to 
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With the approximation Sn = r-- and A„ > \, the partition function reads 
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This result is identical to the result Eq. (11201) to be derived using the spectral density method. 
However, the saddle point derived from Eq. (11061) can be approximated to 



^saddle ~ 2X a ~ 3" 



This approximation in the large N c limit leads to 
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with an additive constant equivalent to the one in the Gaussian-like saddle points method. 
The discrepancy between the Gaussian-like saddle point method and the spectral density of 
the color eigenvalues method comes from a redundant additive constant emerges due to the 
kind of approximation is considered. This redundant constant appears in the gas of only 
adjoint particles in the U(N C ) representation. The SU(N C ) representation does not have this 
kind of problem. In order to verify the results of the Gaussian-like saddle points method, we 
compare them with the results of the method of the spectral density of color eigenvalues Q]. 
The canonical ensemble in the adjoint representation is computed in a similar way done in 
the fundamental representation. Using a similar transformation, the canonical ensemble in 
the term of the spectral density is reduced to 
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where the integral pre-factor C is analogous to that one given in the fundamental represen- 
tation. Using the symmetry property of the spectral density p(9) in the range — n < 9 < tv, 
the stationary equation reads 
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The solution for the range 2A 9 > 1 is already known 
reads 
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However, in the limit X g ^> 1, we get sin 2 (6* /2) = 1/ (2A 9 ). In the context of the spectral 
density method the canonical ensemble, in the large iV c limit, becomes 
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The solution for the energy domain X g < | <C A 9 ^^ max j mum < 1 reads 

\\iZ{X g ) ?a Constant x exp [function (A 5 )] , A 9 < -, (121) 

for the entire color range \9i\ < ir, where the function of X g is independent on N c . Hence, in 
the limit N c — > oo, we have 

lim = 0. (122) 

It means that in the limit N c — > oo, the adjoint particle contribution vanishes and only 
the fundamental particles contribute in the low-lying energy solution for the energy range 
A g < This result is interpreted as follows: the low- lying gluonic spectrum vanishes and the 
glueballs appear only as highly excited states in the high-lying mass spectrum. Fortunately 
the existed hadronic mass spectrum agrees this interpretation. 

The logarithm of partition function lnZ(A ff ) for the gas of color-singlet bags of adjoint 
particles versus the thermal partition parameter X g with various color numbers N c =2, 3, 
4 and 5 is displayed in Fig.flSj). The partition function for the high-lying energy X g > X 
is solved using the spectral density of eigenvalues method and the Gaussian-like saddle 
points method. The both solutions are extrapolated to the low- lying energy range X g < A . 
It is shown that the solution of the Gaussian-like saddle points method is adjacent and 
almost parallel to the density of eigenvalues method and they are splitting by a tiny additive 
constant over the high- lying energy domain X g > |. Furthermore, the exact numerical 
solution is found in the midway between the two solutions. The tiny split between the 
exact numerical solution and the asymptotic high-lying energy solution for the two methods 
becomes less pronounced as the number of colors N c increases. The extrapolation of the 
asymptotic high-lying energy solution which is found by either methods deviates from the 
exact numerical solution over the thermal running parameter runs over the range X g < \. 
This deviation increases significantly as X g decreases far away from the critical point. The 
extrapolation of the high-lying energy solution fails to fit the low-lying energy solution. This 
characteristic behavior indicates a possible phase transition to another analytic solution. The 
high-lying energy solution is simply deflected at the Gross- Witten point. 
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The asymptotic low-lying energy solution versus X g in the energy domain \ g < ^ g r C T\t\cal) 
is displayed in Fig. The asymptotic analytical solution for the low-lying energy states is 
basically obtained by the power expansion over X g and the orthogonal relations in the fun- 
damental group representation. The solution for adjoint representation is obtained by the 
convolution trick and introducing the Lagrange multiplier. The exact numerical solution co- 
incides with the asymptotic low- lying energy solution in the energy domain X g < -^^^critical) ■ 
When the thermal running parameter reaches the Gross- Witten point -^(critical) ~* ^° = I' 
it deviates from the exact numerical one and this deviation grows up significantly as X g 
increases and exceeds the Gross- Witten point Ao = \- The extrapolation of the low-lying 
energy solution fails to fit the high-lying energy solution and subsequently this mechanism 
indicates that the low-lying energy solution is deflected and a phase transition to another 
analytical solution takes place. In order to demonstrate that the analytical low-lying energy 
solution is a correct one, We expand it with respect to the thermal partition parameter 
X g . The results are displayed in Fig.©. It is shown that the asymptotic analytic low-lying 
energy solution coincides with the exact numerical one and when the analytical solution is 
truncated to a lower order it deviates from the exact numerical one significantly. 

VII. FUNDAMENTAL AND ADJOINT REPRESENTATION FOR SU(3) AND 
U(3) 

In the QCD the low-lying hadronic states and the fireballs are treated as bags consisting 
of fundamental quark-antiquark particles and adjoint gluon particles as well. Generally 
speaking, the adjoint particles are assumed to be the interaction among the fundamental 
particles, their anti-particles and the adjoint particles themselves. It will be a reliable 
approximation to treat the hadronic states as an ideal gas of fundamental and adjoint 
particles in the color-singlet state. It would be also a good approximation to ignore the 
interaction among the constituent particles in the present model. 

A. Maxwell-Boltzmann statistics 

For the highly thermal excitations, the Maxwell-Boltzmann statistics becomes a reliable 
one. Fortunately, the Maxwell-Boltzmann statistics is relatively a simple one. In the context 
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of this statistics, the partition function for the fundamental and adjoint particles reads 



Z(\,\q) = J dfi(g)esxp A g ^-tr c R ad j(#) + X qq — tr c (R fuil (ff) + Rf un (#)) 
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In order to simplify our notation we rescale the thermal running parameters X g = , 
dim 9 = Nl — 1, and X qq = With these redefinitions the canonical ensemble is simplified 
to 
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For the large thermal running parameter limit X g > A fl ^ 7 ^ m j n and X qq > \_ ( - // - )m j n , it 
is straightforward to evaluate the partition function using the Gaussian-like saddle points 
method. In the asymptotic large X g and X qq , the saddle points are dominated in a narrow 
interval around the origin. The expansion around the saddle points Q{ ~ reduces the 
partition function to 
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The invariance measure is approximated to 
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The Gaussian-like integration over 5C/ (JV C ) leads to 
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The same calculation can be carried in the same manner for U(N C ) group. In the asymp- 
totic large thermal running parameters X qq > A ? -^ 7 ^ m j n and X g > A g ^^ m j n the canonical 
ensemble for U (N c ) becomes 



Z(ii)(X g , A, 
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where 
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The canonical ensemble is simplified for a further analytical investigation. Using the Stir- 
ling's approximation in the large N c — > oo limit and approximating the summation over n 
to an integration over n, the canonical ensemble can be approximated to 
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In contrary, it is hard to compute the spectral density in the large X g and X qq limits in the 
context of the spectral density of color eigenvalues method. Nonetheless, in that method, 
the canonical ensemble can be calculated only when the spectral density is known at first. 

Under certain conditions when both X g and X qq become small, the color eigenvalues are 
distributed uniformly over the entire interval — 7r < 9i < n. In this case, the partition 
function in the U (N c ) representation is written as follows 
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where the Lagrange multiplier £ have been introduced in order to simplify the above equa- 
tion. By expanding the exponential and evaluating the series term by term using the group 
bases orthogonality and then finally re-summing the resultant terms, the partition function 
is reduced to 
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where £o is determined by the constraint 
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This leads to an additional constraint (q<1. It is possible to approximate Eq. (11331) for the 



energy domain X qq - < | in the following way 
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After evaluating the integral over the parameter £ in Eq. (11351) . the partition function is 
approximated to 
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where 
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The integral on the right hand side of Eq. H136[) can be approximated to 
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Using the approximation given by Eq. (11381) . the partition function is reduced to 
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In order to verify the result, we assume that X g = for a system consisting only of funda- 
mental particles. The partition function is reduced to 
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On the other hand, in the limit X qq - = 0, Eq. (ll39p is converted to 
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for a system consisting only of adjoint particles. 

We are in a position to determine the critical point for the phase transition. This 
can be done by analyzing the canonical partition function. The solution (II) is given 
for the thermal running parameters run over the ranges (\ g > X gCI [t[ ca \ > \(j/)min) anc ^ 
(A* > A 

(^critical — ^qq(ii)min) ■ Usually this solution corresponds the continuous high-lying 
hadronic mass spectrum. The logarithm of solution (II) is given by Eq. (jl3ip . and it reads 
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(142) 



This function is concave up where the minimum is located at the points A 5( - J/ ^ m j n and 
^qq(ii)mm- However, the solution along the thermal running coupling axis in the range less 
than the minimum points X g ^ II - )m ^ n and A 9 -^ 7 ^ m j n is excluded since it will not be physical. 
Usually, the running thermal parameters X qq - and X g are not strictly independent but rather 
they are functions of the variable (. They can be written as A ff (£) and X q q(C)- The variable £ 
could be the hadronic mass m or V/f3 3 etc. The extreme left-side threshold of the solution's 
(II) range is determined by finding the solution's minimum location as follows 
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The minimum thresholds are given by 

\g(//)min = ^99 (C(/7")mm y 
\(//)min = ^9 (C(//)min) ■ (145) 

However, in the case that X qq and X g are independent parameters then the minimum points 
are determined by the set of equations 

l im 9 f lnZ (ii)C\,\g) \ =Q 

lim S( lnZ ^y X ^) = 0. (146) 
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The set of constraints for the solution (II) is given by 

\q — \g(77)min = ^> 
or 

^ 9 — \(//)min = 2' 
(The allowed range for solution (II)) . (147) 

However, the number of the constraints must be kept minimum. In the case the thermal 
running parameters X qq and X g are strictly independent on each other then the the thermal 
running parameter X qq for the fundamental particles is adopted as the master constraint 
for the minimum limit while the constraint for X g is restrained in order to get a feasible 
solution. Nevertheless, in the realistic QCD, the both parameters X qq and X g depend on a 
characteristic variable such as the hadronic mass m or even V/ ft 3 as mentioned above. 

On the other hand, the solution (I) that runs over the range (0 < X qq < Agg(/) ma x) and 
(o < X g < A 9( /)max), is given by Eq.(|139p 

lim lnZ ^y^ „ (148) 
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where \ qq - and X g are again assumed strictly independent on each other. The solution (I), 
usually, corresponds the discrete low-lying hadronic mass spectrum. Furthermore, the set 
of constraints which are associated with the solution (I) reads 

\% < 1, 

A qq ^ A qq(I)m.8iX X ! 

(The allowed range for the solution (I)) . (149) 

The critical points X q q and X g for the phase transition are determined at the point where 
the two asymptotic solutions match each other 

toZ,,,(A,,A«) = lnZ {JI) (X X*) (\ ^) critical ' < 150 > 

Jvf-Large ^ Tvf^Large iv c v /critical 

The both solutions are found match each other at the critical points 

\qcritical ~~ 2' 

\critical = °- ( 151 ) 
This solution is very interesting and is interpreted that the low-lying hadronic mass spectrum 
does not mix the fundamental particles and the adjoint ones significantly. The hadronic 
states with the constituent fundamental particles (e.g. quarks) has no significant adjoint 
components in the low energy domain. This scenario explains why the low-lying hadronic 
mass spectrum is likely consisting of mesons and baryons and is unlikely consisting of exotic 
hadronic states mixed with significant gluonic components (e.g. hybrid: quarks, antiquarks 
and glueballs are not common in the low-lying mass spectrum). The second example is the 
QCD with N c = 3 ^fl avor = 2. According to Eqs. (11571) and fl!64j) . we have \ g = ^\qq- In 
this example, we get from the constraint given by Eq. (11501) the critical point -^gg CI -itical ~ 
0-6 > A 9 -/ Jr x m j n where ^^f/nmin ~ 0-18. It is shown that the low-lying solution given by 
Eq. (1148p has not an adjoint component and subsequently it leads to an over estimation for 
\qcriticah H° wever f° r a rather large adjoint component the solution given by Eq. (11481) is 
not appropriate to determine precise the critical point. In order to determine the threshold 
limit of the (pre-)critical point we can return to Eq. (11251) . It is noticed that Eq. (1125j) can 
be written as a semi-fundamental representation as follows 



4V)/^K« 2 ^> E "' 9 ". (152) 
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The (pre-)critical point is determined by 

A^^(2JV A 55 + 2JV c A 9 ) cr . t . cal = i. (153) 

This result is analogous the critical point derived in section IVIH The advantage of this 
method, the explicit expression for the low-lying solution is not needed and it gives the 
threshold limit of the critical point ^(p re _)critical — -^critical - Hereinafter, we shall adopt 
the above procedure to estimate the location of the critical point and set the correspondence 

^ (pre-) critical — > ^critical- 

The grand canonical potential for the color-singlet bags of fundamental and adjoint par- 
ticles versus the thermal running parameter A with various color numbers N c = 2, 3, 4 and 5 
is displayed in Fig.®. In this toy calculations, we have set A = X q g = X g . The asymptotic 
analytical high-lying energy solution is displayed for the Gaussian-like saddle points method. 
The asymptotic high-lying energy solution is compared with the exact numerical solution. 
It is shown that the asymptotic analytical solution fits precisely the exact numerical one for 
the energy domain A > A cr j^j ca j even with a small number of colors N c = 2. The extrapo- 
lation of the high-lying energy solution to the range A < -^ cr itical deviates from the exact 
numerical solution. This deviation becomes significant as A decreases. It is evident that the 
extrapolation of the high-lying energy solution to the low-lying solution is unphysical. More- 
over, it is evident that the grand potential is deflected at the Gross- Witten point. Hence, 
the high-lying energy solution undergoes a phase transition to another analytical function 
that describes the low-lying energy solution more successfully. 

The asymptotic analytical low-lying energy solution is displayed in Fig. (Q. The low-lying 
energy solution is approximated by expanding the canonical ensemble as a power expansion 
over A and then integrating it over the entire color circle and using the Lagrange multiplier 
and the convolution for the adjoint particles. The low-lying energy solution fits the exact 
numerical one over the range A < A cr ^ ca j even with a small number of colors N c = 2. It 
deviates from the exact numerical solution when A increases and exceeds the critical point 
\q+\ g = 2A > 2A cr j|.j ca j. The numerical calculations shows that the split between the exact 
numerical solution and the low- lying energy solution (i.e. solution (I)) becomes noticeable 
at the point A = ~ 0.34 > ^(p re _) cr iti C al = I * n agreement with the result of Eq. (11501) . 
Furthermore, the extrapolation of the low-lying energy solution to a larger thermal running 
parameter A > A cr j^ ca ] deviates from the exact numerical solution significantly. This failure 
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for the low energy extrapolation indicates that the analytic solution is deflected in the middle 
between the low-lying and high-lying energy solutions and proves the existence of the phase 
transition. 

In order to demonstrate the phase transition and the deflection point for the analytic 
solution, we display the low and high-lying energy solutions and their extrapolations in 
Fig. (fTUl) . The asymptotic analytical solutions and their extrapolations for the entire energy 
domain are compared with the exact numerical one. It is shown that the low-lying energy 
solution fits the numerical one for the small thermal running parameter A < --^critical ' 
When the thermal running parameter A reaches the Gross- Witten critical point it starts to 
deviate from the exact numerical one significantly as A increases and runs far away from that 
point. On the other hand, the high- lying energy solution fits the exact numerical one for 
the large thermal running parameter over the range A > A cr j|-j ca j while its extrapolation to 
small values deviates from the exact numerical one and this deviation increases significantly 
as the thermal running parameter A decreases down below the Gross-Witten point and 
approaches the origin. Evidently, the extrapolations of both solutions fail to fit the exact 
numerical one when they go far away beyond the critical point in the opposite directions. 
This demonstrates how the analytical solution is smoothly deflected at the critical point. 
This soft deflection causes the asymptotic analytic low and high- lying solutions to lie just 
above each other and reduces the intersection possibility to the minimum in any realistic 
rough approximations. Therefore, Gross-Witten point is then determined in the midway 
between the smooth interpolation from the low to high-lying energy domain or from a 
small to the large thermal running parameter A. The exact numerical solution matches the 
asymptotic low-lying energy solution from the left and then it matches the asymptotic high- 
lying energy solution from the right beneath this small interval around the Gross-Witten 
point. However, in the case that the low- lying energy solution is not known, the extreme left 
hand critical point is determined by the minimum point of the high-lying energy solution 
extrapolation. This point is the extreme left hand threshold of the high-lying energy solution. 
This threshold point can shed the light on the existence of the phase transition even for a 
complicated realistic physical situation can not be solved using the spectral density method. 
The characteristic solutions of the low-lying and high-lying energy domains are found not 
restricted to large N c — > oo limit and is satisfactory even for a small number of colors N c = 2. 
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B. Quark and gluon fireball 



The hadronic states are bound states of quarks and gluons confined by a color-singlet 
state. The excited hadronic states are assumed to be Hagedorn states. The Hagedorn states 
are approximated as bags of ideal quark and gluon gases and each bag is in the color-singlet 
state. In the realistic QCD, the quarks are fundamental particles and are satisfying the 
Fermi-Dirac statistics while the gluons are adjoint particles and are satisfying the Bose- 
Einstein statistics. The canonical potential in the Fock space for quark and anti-quark with 
an internal color structure is formalized as follows 



\nZ qq ((3,V) = -Ltr c £ (j> [l + R(</)e"^] + J> [l + R*(g)e 



N 

c Flavors v a 

= w tTc E E ln I 1 + ( R Q?) + R *Q?)) ^ Eqa + z~ WEqa ] ■ ( 154 ) 

c Flavors a 

For the highly thermal excited hadronic states, the color saddle points become dominant 
over a tiny interval around the origin. This highly thermal excited states correspond the 
energy domain A > -^critical — -\//)minimum an< ^ the ana l°gue solution Z(j/)(/3, V). In the 
Hagedorn Hilbert space, the quark-antiquark canonical potential is approximated to 

lnZ q - q (P,V) « ]T ^ln[l + 2e-^+e- 2 ^" 
Flavors a 



1 1 ^ c 
- — V# 2 

c 7. 



= ^-^S ) EC (155) 



Flavors a 1 1 

N c 



where the following parameters are introduced in order to simplify our notations 

^ = E E ln I 1 + 2e ~ l3Eqa + e ~ WEqa > 

Flavors a 

( 2 ) 1 2e~ l3Eqa 

= W CT? ^ ? fl + 2e-^ + e- 2 ^-] ' (156) 
Flavors a 1 1 

The sum over states for an ideal gas of massless quarks and antiquarks confined in a cavity 
with a sharp boundary is approximated by the integration over the density of states as 
follows 



*g> = (2j + l)-iV c n flavor / —2 In 



^(2j + l)iV c n flavor ^, (157) 
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and 



1 r Vd 3 p 2e~ l3p 



^^(2 J + l).iV c n flavor -/_| 



(2tt) 3 (1 + e -ft>) 2 



V 



-(2j + l)n flavor — , 



(158) 



where V is the bag's volume and (2j + 1) comes from the spin degeneracy. Since the 
constituent adjoint gluons are satisfying Bose-Einstein statistics, the canonical potential for 
the gas of gluons reads 



\nZ g ((3,V) = -^tr c £ln[l-R adj (s)e-^- 



(159) 



'-g a 

The internal color structure for SU(N C ) (or U(N C )) symmetry group is separated and eval- 
uated as follows 



oo 

tr c {ln[l-AR adj (s)]} = ]T A*tr c R adj ((? fc ), 



k=i 

oo 

k=i 

oo 
k=l 



cos k {6i — 6j) 



E 



£(X> [l-Ae^" e 



+ (N C -1)}. 
+ (iV c -l)J, 
+ (N c - 1) ln[l - A]. (160) 



The approximation of the gluon canonical potential around the color saddle points, which 
are dominated near the origin, becomes 



N c N c 



In Z g ((3, V) « -^ln[l- e-^-] - ^ £ £ (*„ - * 



n=l m=l 
2 



-(3Egc 



1 ^ e 
dim fl V [1 - e-^° 



^ 0) - E E (fn - e m ) 



(161) 



n=l m=l 



The following terms 



and 



^ 0) = EMl-e- 



9 A\™ ^ 



(162) 



dim s a [1 - e~^ E ^\ 



2 ' 



(163) 
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are introduced for our convenience in order to simplify our notations. Furthermore, we have 
the dimension dim s = (iV c 2 — 1) and the color constraint 9 Nc = — X^f 1 ^ f° r the SU(N C ) 
group representation. Since the quarks and gluons are assumed to be confined in a cavity 
with sharp boundary then the summation over the energy states is performed as follows 



$f = (2j + l)dim 9 



7T 

90 



Vd 3 k 
V 



In 



1 



and 



$( 2 ) 

9 



(2j + l)dim^ 3 , 



(2 ^ + 1) /?27 



(164) 



(2tt) 3 [1 - e -/»] 



2 > 



1 



V 



(165) 



The partition function for an ideal quark and gluon gas, which is in the color-singlet state in 
the Hilbert space of Hagedorn state, is the tensor product of Fock spaces for the quark and 
gluon partition functions. The resultant tensor product is then projected into a color-singlet 
state as follows 



Z q ~ qg (P,V) = J dpL(g)Z g (P,V)Z q - q (P : V). 



(166) 



The integration over the Fourier color variables is evaluated by the Gaussian quadrature 
integration around the saddle points. In the extreme condition in the limit (3 — > 0, the 
color saddle points become dominant in a tiny range around the origin. Hence, the canon- 
ical ensemble for the highly thermal excited states is then approximated by the quadratic 
expansion around the saddle points as follows 



z g - gg (P,v) 



J qqg{U) 



Z q qg(II)(P, V), 

exp [*£> + $r] 

x I dfi(g) saddle ex V 



.I$( 2 _) w_ 

2 qq V 1 2 2N, 



l(2NM)) n c n c 



n=l m=l 



(167) 



where the approximate Haar measure near the saddle points is embedded in the Gaussian 
integration. After evaluating the Gaussian-like integral over the Fourier color variables, the 
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canonical ensemble becomes 

Zqqg(H)(P,V) = Mggg 



2iV c $l 2) + $ (2) 



exp 



$W + $(o) 
m ~ g 



qq 



(168) 



where the pre-factor normalization is given by 

\f 1 (27r) jVe/2 nf=i fcj fl69 . 

In order to be in a position to compute the micro-canonical ensemble, it is appropriate to 
extract the thermodynamical ensembles (V, (3 3 ) in the following way, 

s /V 



Z qqg{II) ((3,V) = N q 



/Q3 V 9« 9 / 







(170) 



The functions ${9, $(°), and $( 2 ) are independent on V and (5 and they are calculated 
by dividing the un-tilde terms by J^, 



qq 




qq 


qq/ 03 


$(°) 

9 


9 / p 3 


$( 2 ) 
S 


9 / /3 3 



(171) 

The micro-canonical ensemble is found by calculating the inverse Laplace transformation 
as follows, 



Z q - qg{II) (W,V) 



2ni Jpc-ioo 



d^ w Z qqg{u) {^V), 



qqg 



2ni Jpc-ioc 



(&°J + $(o) 
V qq ~ g 



jv-?- 



T- ( 172 ) 



It is appropriate to define the following constants explicitly which are independent on (5 and 
V, 



(o) 
qqg 



(°) j_ <f>(0) 
7tt 2 



qq ~ 9 



= (2j + 1) 



7T 



360 



^flavor + 7^(^-1) 



(173) 
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and 



^+2iV c $f = ^(2j + l)U avor + iV ( 



6 



174) 



The Laplace transform is evaluated using the steepest descent method. The Laplace saddle 
point is determined at the following point 



Km' 

w qgg 



(175) 



The micro-canonical ensemble for the Hagedorn states corresponds the high-lying energy 
solution as follows 

%agedorn(^) ~ Z q , g(II) (W,V). (176) 
The micro-canonical ensemble for quark and gluon bag in the color-singlet state reads 



J qqg{H) 



(W,V) 



2>/2F q " *W \VW 3 



cxp 



;i77) 



where the pre-factor constant is given by 

V ma) 



(2NM ] + ) 



N2_ 



(178) 



After scaling the bag energy with respect to its volume x = W/V which is denoted as the 
bag energy density, the micro-canonical ensemble becomes 



Zqqg(II)(x,V) = 



1 



cxp 



(3*&)V/V 



(179) 



2V27T 

In the real world of the relativistic heavy ion collisions, the story is more complicated. 
The masses for the constituent quarks even for the light flavors do not identically vanish 
to zero under the Gross- Witten phase transition line to the Hagedorn phase due the chiral 
bound interaction. Although the constituent zero mass approximation for the light fla- 
vors is satisfactory in particular for the hadronic phase just below the deconfinement phase 
transition point, the constituent strange flavor mass still deviates significantly from zero. 
Therefore, the inclusion of massive flavors such as strangeness degree of freedom will dras- 
tically modify the numerical results quantitatively. The calculation of the micro-canonical 
ensemble becomes more complicated for the massive flavors. The major trouble comes from 
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the analytical calculation for the location of the Laplace saddle point in the steepest de- 
scent method. In this case, after the integration by parts the quark and antiquark partition 
function namely, 



*SU*S(AV), 



*£> W, V) = (2j + l)iV c n flavor / ^2 In [l + e "*^ 



(180) 



is reduced to 



*<g(P,V) = (2j + l)N c n & ^ 



VI r°° „ 1 /W 



VOT (5* 7T 2 7o 3 #FT^? 



V 1 



(2j + l)^fla V or^^V.»W)' 



;i8i) 



where the integration term reads 



. (o) (m/3) = / 

99 JO 



oo I 

dx 



X 



3 \/x' 2 + (5 2 m 2 



(182) 



The same thing can also be done for the second order function When the mixed- 



grand canonical ensemble is transformed to the micro-canonical ensemble, the massive flavor 
produces a complicated transcendental function for the stationary Laplace saddle point 
solution. Generally speaking, it is possible to solve the transcendental equation by the 



iteration. At first, the solution for the zero^ n iteration is assumed to be (5^'m = 0. Then 
in the first iteration, we solve the stationary saddle point (3q = (5^p . In the second iteration, 
the massive flavor is considered explicitly and then we repeat the iteration to the higher 
order until the convergence is achieved. This procedure is summarized as follows, 



(°)„ 



(5m 
(5m 
(5m 



0, 

(5ipm, 



(183) 



Nonetheless, it is expected that the convergence is achieved very rapidly and subsequently 
the first order iteration is an adequate and sufficient approximation for light flavors including 
the strangeness. Therefore, the first order iteration truncation (5m — > (5^m is an appropriate 
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one for the relativistic heavy ion collisions 

V) = (2j + l)7V c7lflavor ^l/ $(0) (/#> m ) . (184) 

The value of the point /3® can be simplified in the context of the standard MIT bag model. 
After considering the bag's total energy m = W + BV and the volume V = m/AB con- 
straints, the stationary point becomes 



= ^- (185) 

The term Tp can be interpreted as an effective temperature for the screening mass. Alter- 
natively, in order to overcome the possible complication in the micro-canonical ensemble, it 
is possible to assume the following ansatz 

Tfl 

(3m -> — , and T(3 ~ 0(1), (186) 

for the high temperature approximation. This assumption simplifies the mixed-grand canon- 
ical ensemble significantly, 

V) = (2j + l)iV c n flavor ^/ $ („) . (187) 

An alternate assumption regarding the highly compressed hadronic matter reads 

(3m oc -m. (188) 
[i 

Nonetheless, the density of states for the highly compressed matter at low temperature needs 
more consideration. In this regime the physics is rich due to the configuration space, flavor 
and/or flavor-color correlations. Furthermore, a new equation of state is expected to emerge 
due to the formation of the color superconductivity. 

On the other hand in the limit of low-lying energy (3 1 (i.e. the diluted and relatively 
cold matter which is analogous to Z^(X) solution where A < -^critical — ^OOmax), the 
partition function for an ideal gas of quarks and antiquaries can be approximated to 

lnZ q - q (P,V) = tr c £ £m[l + R(2)e-^-]+tr c £ £ln [l + R*(^)e"^ 
flavor a flavor a 



E 

flavor 



tr c R(0) ^ + tr c R* (g) ]T 

a a 

E (Ee-^)tr c [R(^) + R*(^)]. (189) 



flavor 
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The low-lying energy limit for the gluon gas reads 

\nZ a {p,V) = -tr^lnfl-Radj^e-^-], 
« ^E^ad^) 6 ^ 

a 

= (E^ Q )tr c R adj (^)- (190) 
By substituting the following thermal running parameters 

\- q = E E^-, 

flavor a 

K = E e "^ 9Q ' ( 191 ) 

a 

we get the partition function 

WA* = f *( 9 )/* a! ^ t^ • R(S " ^ • R ' <^ V"* tr • ,R(9,+Il • <9,, , (192) 

J— 7T 

for the SU(N C ) group representation. This equation is similar to Eq. (1132p . The result reads 

Z q q g (i)(X q q, X g ) = - 1 Ag exp [\%/ , (193) 

where X q g = ^ ^^^~ n flavor ^ or mass l ess flavors and A 9 = ^ Nc n 2^ ^ 2j '^ y • This approxi- 
mation is an appropriate solution over the energy domain X qq - < -^g^critical — \g(/)max 
and X g < -^^critical — \(-0max- This result is analogous to the result given in Eq. (11391) . 
Furthermore, we have for massive flavors the following thermal gg-running parameter 

(2i 4- 1)V m 2 B 2 

\« = ^flavor ^3 ^ Mm q i3), (194) 

where K 2 (x) is Bessel function of second kind. In the case of the low-lying energy with 
massive flavors in the limit (3m q ^> 1, the gg-running parameter becomes 

\,~^jM±iK^y\-^. (195) 

Nonetheless, in the realistic case, the density of states derived from the micro-canonical 
transformation of the low-lying mixed-canonical ensemble given by Eq. (ll93p will be found 
in Section IVIIII to be replaced by the discrete low- lying mass spectrum. 
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VIII. MICRO-CANONICAL ENSEMBLE AS A DENSITY OF STATES 



The micro-canonical ensemble becomes known after calculating the inverse Laplace trans- 
form of the mixed-canonical ensemble, 



Z(W,V) = — d^ w Z(f3,V), 

ZlTl J/3 c ~ioo 
1 rl3 c +i<x> 

= — dp exp\J3W + hi Z(P,V)]. (196) 

2m J/3 c -ioo 

The above integral can be computed by the steepest descent method. In this method the 
integral is approximated to the Gaussian-like integral around the saddle point. The general 
result reads 



1 1 exp 

Z(W,V) = ' 



(3W + In Z(p, V) 



2^ 



( _ l dZ(/3,V) \ z _ d 2 Z(/3V) ] 

\Z{(3,V) d(3 ) ZQ3.V) W 



(197) 



(3=/3 

where the saddle point (3 of the stationary condition is determined at the extremum point 



= 0. (198) 



Z(f3,V) dp 

In the standard MIT bag model, the density of states is determined for constituent particles 
confined in a cavity with a sharp boundary. In a realistic model, it is possible to deform 
the cavity's boundary. The (^-function for the sharp boundary in the standard MIT bag 
model can be smeared by the Gaussian smoothing function. The extreme conditions in the 
relativistic heavy ion collisions may smooth the sharp boundary for the quark and gluon 
bag. Therefore, the micro-canonical ensemble for the gas of bags with extended boundaries 
is approximated to 

Z(W,V) = J dv8(v-V)Z(W,v), 

= J ^/smooth (v-V)Z(W,v), (199) 

where the boundary surface is extended by the Gaussian smoothing function 

; e -#^-^ 2 . (200) 



^smooth — ^ 



A/W 

The bag's boundary becomes more extended with respect to the energy when the bag's 
energy increases with respect to temperature, while on the contrary, the surface turns to be 
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sharper for the low energy. The bag volume is proportional to its energy vq oc W. In the 
low energy limit, the bag's boundary is reduced to a 5-function 

|™ /smooth^) = S ( v ~ v o)- ( 201 ) 

Hence, the Hagedorn's density of states for the standard MIT bag with a sharp surface 
boundary reads 

(in \ 
v-—)Z(W,v), (202) 

where W = m — Bv. The micro-canonical ensemble Z(W,v) is given by Eq. ( 1177ft . The 
density of states can be re-written as follows 

p {II) {m,v) = 5{m-ABv)Z(m,v). (203) 

The generalization of the micro-canonical ensemble to take into account the inclusion of 
volume variation reads 

P(ii)(W,v) = f smooth (v)Z(W,v). (204) 

The volume fluctuation effect is studied in details in Ref.j28] and the bag stability needs a 
further investigation. 

A. The critical mass for the bag consisting of fundamental particles obeying 
Maxwell-Boltzmann statistics 

The canonical ensemble for a gas with the internal color symmetry U(N C ) reads 

Z((3,V) « Z (I) (P,V), 
Z (I) (J3, V) = e N *~ x \\ < l -< A (/)max « 1, (205) 

and 



Z(P,V) « z {II) (p,v), 

Z m (f), V) = (i) 2 eK'M".') (6 A > \ > A„ J)min ) , (206) 
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for the low- and high-lying energy solutions, respectively. The effective thermal running 
parameter A for an ideal gas of fundamental particles embedded in the thermal bath and 
satisfying Maxwell-Boltzmann statistics reads 



A = A/ATf, 

_ .1(27 + 1) [Yfp e -PV^ 

1 (2j + 1)V 



m=0 



(207) 



7T 2 N c /3 3 ' 

The low temperature 1/(3 phase has a free energy of order one 0(N®). The high temperature 
phase has a free energy of order O(N^). It is characterized by the Hagedorn growth in its 
density of states. 

The micro-canonical ensemble is computed by calculating the inverse Laplace transform 
as follows 



Z {H) (W,V) = lim — / dPe^ w Z (II) (P,V) 

W^oc 2m Jf3 c -ioo 



— 2~ (^) m i n 



2V27T V12iV, 



x exp 



n 2 N r 



^3iV c 2 +5)/2 v ,_(iv2 + i)/2 
1/4 



2(2j + l)J 
4/ 6(2i + l)7V c \v- _3 
3\ n 2 J 4 c 



(208) 



The stationary point for the steepest descent method is found at 

fQ(2j + l)N c V\* 



7T Z 



w 



(209) 



In the micro-canonical representation, the effective thermal running parameter and the en- 
ergy constraint for the highly excited states are reduced to the following 

1/4 



— 2 — \^)min- 



(210) 



In the standard MIT bag model, the effective bag's mass and volume are given by the 
relations 



m = W + BV, 



(211) 
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m 
W 
V 



ABV, 
3 

— m, 
4 ' 
m 

4B' 



(212) 
(213) 
(214) 



In the present model, the Hagedorn states appear in the hadronic mass spectrum when 
the energy reaches the threshold A > (pre-) critical — 2 — ^(//)min an< ^ subsequently the 
Hagedorn mass threshold exceeds the limit 



m > 2v / 2TiVi 1+3/4) 



2B 



.1/4 



,(2J + 1) 

> 6.2GeV with B 1/A = 180MeV. 



(215) 



This threshold mass is relatively large due to the classical Maxwell-Boltzmann statistics that 
is considered in the present scenario. 

On the other hand, the situation is rather different for the low-lying energy for the small 
thermal running parameter A < \ < A(/)max (usually A(/)max < 1), where the micro- 
canonical ensemble is evaluated as follows 

— / d[3e? w Z {I) ([3,V), 

c+ioo \,pvl. 

(216) 



Z(i)(W,V) 



— f d(3e? w e 
2m J/3 c -ioo 



where d„ 



(23+1) 



. In the standard MIT bag, we have the effective bag energy W 



Tin. 



For the energy excitation less than the threshold A < \ < A(j)max> we have the low-lying 
hadronic mass spectrum and this low-lying spectrum ends when the energy scale reaches the 
critical point 

: 6.2GeV, with B 1/4 



m < m A 



(pre-) critical 
< m(/)max : lOGeV, 



180MeV, 



(217) 



where the upper limit m^max will be given by Eq. (12511) . When the energy scale exceeds 
this point, the micro-canonical ensemble for the low- lying hadronic states becomes an inap- 
propriate. The extrapolation of the steepest descent approximation leads to the density of 
states, 

1/14 



Z (I) (W,V) 



1 



6dl 



2 v / 2^v / 7V r ~ 1 /7H M / 7 



cxp 



V^W-^exp [(• ■ .)V 2 "W" 7 ] (m < m (/)max ) . 



(218) 
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This approximation is not valid for small W or m < mmmax if the low-lying spectrum is a 
discrete one as will be found in the realistic physical situation in section IVIH B[ We shall 
assume that m(/) max > "^Critical ~ m physical w h ere m physical ~ ^ ® GeV is the maximum 



mass for the known physical particles found experimentally without strangeness [45J. For the 
low-lying hadronic mass spectrum, we have considered the density of states for the known 
mass spectrum particles 

^Low-LyingW V ) " Z H)(W,V), 

PLow-LyingW^) <* ^Low-LyingW y )> 

m Critical 

ps ^ 5(m — mi)5(v — Vi). (219) 

rrii,vi 

Nonetheless, the micro-canonical extrapolation given in Eq. (12181) is useful to show the 
discrete hadronic low-lying mass spectrum is essential. If the standard MIT bag model 
with a sharp surface boundary is applied to Eq. (|218p . no phase transition to an explo- 
sive quark-gluon plasma takes place. To see this let us adopt the standard MIT bag ap- 
proximation where W oc V. This consideration leads to a continuous density of states 
p ~ (• • ■)v~ 3 / 7 e("^ v8/7 . This density leads to a divergence in the isobaric partition function 
and this standard approximation of the MIT bag for the low-lying mass spectrum fails. It 
is evident that the Hagedorn states do not appear in the low- lying mass spectrum limit. 
Furthermore, if we argue that the continuous density of states grows exponentially as a 
linear exponential growth with respect to the hadron's excluded volume, it is possible to 
extrapolate V 2 ^ 7 W 6 ^ 7 oc V or W oc V 5 ^ 6 . With this approximation, the density of states is 
reduced roughly speaking to p ~ (• • ■)v^ 1 ^ 3 e^'"^ v . This density gives no phase transition to 
the quark-gluon plasma. The both above approximations show that continuous low-lying 
mass spectrum can not generate explosive bags. 



B. The critical mass for the quark and gluon fireball (e.g. Hagedorn state thresh- 
old) 

In QCD, Hagedorn states are assumed to be a bag of weakly interacting quarks and gluons 
confined in a color-singlet state. The canonical ensemble for a quark and gluon gas in the 
color-singlet state is given by Eq. (1170j) . The logarithm of the canonical partition function is 
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simplified as follows 

log [z qqg (v/(3 3 )) « log (z qqg{II) (v/(3 3 )) (for the energy domain X qqg > A 



qqg(II)m.mj 



\og(z mg{n) (v/p)) = 1 



/3 3 



99 ^ 9 



iV c 2 -l, 7 
— - — log — + constant, 
I p 6 



(220) 



for the energy domain X qqg > \q g ni)mm wnere ^999 ^ V/P 3 - The critical point for the 
phase transition is determined by finding the extremum left-hand point as follows 

d 



logZ q - qg{II) (v/[3 3 )=0. 



(221) 



The parameter (j^J is related to the running coupling parameter A and it determines the 
threshold point of a possible phase transition to the Hagedorn phase. The extreme left-hand 
side location for the prospective phase transition is located at 



(0) 

99 



$(°) 
9 



> 



(222) 



Below the minimum left-hand side point, the solution must be already deflected to another 
analytical solution and in this sense the critical point is bounded from below by the extreme 
left-hand side point. The minimum limit for the threshold Hagedorn mass (i.e. ^^(//jixiiu 
where ^ cr iti ca l > m (//)min) ^ or ^ ne P nase transition from the low-lying mass spectrum to 
the highly excited spectrum of fireballs or more precisely the continuous Hagedorn states in 
the context of the micro-canonical ensemble is given by the condition 



3 (<£$ + $f) vw* 



1 

> -. 
- 2 



(223) 



3(iV2 _ i; 

The stationary point of the steepest descent method for the inverse Laplace transform to 
the micro-canonical ensemble is determined by Eq. (11751) . In the context of standard MIT 
bag model with a sharp boundary surface, the point of the critical point is determined by 
the constraint 

1/4 3m 



3VW" 



(224) 



The minimum limit for the threshold Hagedorn mass production in the context of the MIT 
bag with the sharp surface boundary approximation is given by the following constraint 



1/1 



7n 2 7T 2 

N e na awxr H (iVA 2 - 1) 

360 tiavor 90 V c ; 



1/4 



gj + 1) 
2(iV c 2 - 

— >• m > ^(/n m in (the resultant constraint). 



m 



> 1, 



(225) 
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Hence, the threshold Hagedorn mass production for N c = 3 and B 1 / 4 = 180MeV is given by 



m 



> 2206MeV(one massless flavor), 

> 2029MeV(two massless flavors), 

> 1906MeV(three massless flavors). (226) 

The critical point for the phase transition from the discrete low-lying mass spectrum to the 
continuous Hagedorn density of states is located just above the maximum mass of the known 



hadronic particles found in the Data book[45j]. The highest experimental non-strange hadron 
state is roughly estimated to m max ~ 2.0 — 2.3 GeV. The critical point exists above the 
minimum point solution of the Hagedorn threshold ^ cr iti ca l — m (ii)mm- The estimation 
of the Hagedorn's critical mass threshold will be given in Eqs. fl270p and (1271 p . 

The low-lying mass spectrum is considered a discrete density of states that includes all 
the known mass spectrum particles 



^Low-lying(^) = Z^ g{I) (W,V), 

m max 

^Low-lying ( m ' u ) ~ S ( m ~ m i) 5 ( v ~ v i)- ( 227 ) 

rrii,Vi 

Using the following connections between the canonical ensemble and the density of the states 

p(W,V) oc Z qMI) (W,V) + Z (qMII) (W,V), 

« ^Low-LyingW V) + % ire balls(^ V), (228) 

the density of states for the entire hadronic states is approximated to 

m max 

p(m,v) ~ 5{m — mi)5{y — Vi) 

rrti,Vi 

+ ^Fireballs v ) x ( vo lume-mass fluctuation), (229) 

where the micro-canonical ensemble Z q g g nn(W, V) is the continuous mass spectrum for the 
Hagedorn states (i.e. -^Fireballs ^))' ^ n ^ ne standard MIT bag model with the sharp 
surface boundary the volume-mass fluctuation is determined by the relation S(m — 4BV). 
For any fuzzy bag model, the volume-mass fluctuation becomes nontrivial and may change 
the order of phase transition. 
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IX. DENSITY OF STATES 



The density of states can be evaluated by using the phase-space integral. The calculation 
for an ideal gas of particles in finite size has been evaluated and studied without the internal 
color structure {3} . We use the same procedure to re-derive the density of states reached in 
the previous sections in particular section IVIIBI The density of states from the phase-space 
integral reads, 



iV c n (2j + 1)" 



n=2 



rv. 



n 

i=l 



Vd 3 Pi 



(230) 



i=l 



where the constituent particle's energy e(pi) = ypf + m? and it is reduced to e(pi) = \p\ 
for the massless one. Here j is the particle's quantum number and V and W are the bag's 
volume and energy, respectively. The internal color symmetry is embedded in the phase- 
space integral as follows, 



p(w,v) = £ 



(2j + l) r 



n=2 

oo 



ntr[R(p) + R*(p)] 



li=l 



Vd 3 Pi 



i=i 



n=2 



nl 



1 1 

= W^ 2 n\(3n-iyi a{W3V}{K) 
In Eq. (123 ip . the following function has been introduced 



a {W 3 V} (R) = a {w : iv} tr[R(g) + R*(g)} 
i {W 3 V} tr ( R(g) J 



(231) 



n 2 ds 



2tt 2 



where 



a {W :i V} 



(2j + 1) 



7T Z 



x 2 e~ x 



W A V 



dx 
2^2 



x 2 e~ x 



(232) 



(233) 



The summation over the number n is approximated with the help of the Stirling's formula 
to an integration over the variable n 



1 



p(W,V) = — J dnexp 



n J 



(234) 



In (a^w3yy(Ei)\ — Inn! — ln(3n — 1)! . 

The above integration is evaluated using the saddle point approximation, the saddle point 
is found as follows 

i 

a-{H/ 3 y}(R') \ 4 



n 



27 



(235) 
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In this approximation, the density of states is reduced to 



256 /X3 xl/4 
— a {w s v} (K 



xC n {W,V), 



(236) 



where 



C n (W,V) 



o 



dnexp 



Ir— ( a {W 3 V} (R 
2 V 256 x 27 



1 / 256 x 27 

2 \a{ W /3 V /}(R) 

1/8 



1/4 



n 



(237) 



Hence the density of states reads 

1 



p(W,V) = —j= {3a {W 3 V} (R)Y ^exp - (3a {W 3 V} (R 



W 



1/4' 



(238) 



After tracing over the color index for U(N C ) (or SU(N C )), Eq. f)232p becomes 



fl{iy 3 y}(R-) = 2 01^3 cos 9j 



(239) 



i=l 



The projection over the color-singlet state reduces the density of states to 



/'Singlet 



dp(g)p(W). 



(240) 



The density of states under the integral is approximated to the Gaussian-like integral over 
the color variables since it dominates a tiny range around the origin along the real axis of 
Fourier color variables 

\ 1/41 



p(W,V) 



1 



iV:: 



4V2tt 
C p x exp 



6a{ W 3 V } ^2 cos Oi 



w exp 



4 / n c 
- \ 6a {W 3 V} J2 cos 9i 



~2 V ^ 

i=i 



(241) 



where the pre-exponential coefficient reads 



1 



GN r a 



4V27T 

The exponential term is given by 

V„ = 



c Ct{W 3 V} 



W 



exp 



(QN c a {W 3 V} 



i/4- 



(242) 



4 1 / R - \V4 

i i 6^+m w s v 



3N r 



.1/4 



7T 



(243) 
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The color-singlet density of states for the groups U(N C ) and SU(N C ) is, respectively, approx- 
imated to 

PSinglet^^) = C p J ^ ex P -2 V Pl2 e 

c L i=l 

N*/2 

forC/(A c; . 

(V c 2 -l)/2 

! x C p , for SU(N e 



n^n! 



iV c !(2^/ 2 X ^J XCP ' iOTU(Nc 



n^n! /2tt 



N c \(2tt) n °/ 2 )I N c 
(@ V p /N c >\). 



V, 



(244) 



Furthermore, the approximation which is carried out in Eqs. (12411 ) and (1244ft requires the 
constraint T> p /N c > 1. This is derived trivially using the analogous behavior p(A) oc 
e 2N -~ x J d^(g) exp [-±(2iV c A) Ei ^1 where A > 1/2. The constraint V p /N c > 1 is exactly 
the same constraint given by Eq. (l210l) . The asymptotic color-singlet density of states for 
the U(N C ) and SU(N C ) group representations behaves, respectively, as follows 



PSinglet^^) K 





■ 1 ■ 


c 
8 




[w 3 v. 


exp 






iv2-2 




- 1 ■ 


8 




[WW. 


exp 


(@ 


V p /N c > 1) . 



constant 



constant 



(wV) 
(w 3 v) 



1/4- 



1/4 



for Z7(JV C ), 
for SU(N C ), 



(245) 



In order to analyze the results for the low energy limit, it would be worthy to do the 
analysis in the language of the canonical-like ensemble rather than the micro-canonical 
ensemble. The density of states in term of the canonical-like ensemble can be defined as 
follows 



^Singlet ( W > V ) = ^^^( 3a ^ 3y > (R) ) 8exp ^( 3a {^ 3 n( R ^' ' 



1 1 1 



d^e^ exp 



c 



2 W 2m J( - 

The projection of the color-singlet state reduces the density of states to 

1 1 r 1 /-Co+*oo 



(246) 



Jy=l 



^ ( a {w 3 v}(R-" 

/J 
dn(g)exp — (a {w s v} (R] 
. C 



1 1 1 



Co+ioo 



j/=l 



y=l 



(247) 
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Hence, the analysis of the above equation becomes similar to the Gross- Witten critical point 
solution. In the case <C 1 the color saddle points will be distributed uniformly over 

the invariance measure J^ n d/i(g). It is evaluated as follows 



PLow-Lying^ V ) 



111 Ko+ioO f-TT 

__ / dCe Cy 

2 W 2ni Jco-ioo 



dfx(g) exp 



1 {W 3 V} 

C 3 



(trR(p) + trR*(^)) 



y=l 



1 1 1 



dC,e^ y exp 



a{w 3 v} 

, c 3 , 



y=l 



2 W 2iri . 

(i.e.: continuous low- lying spectrum). (248) 

Under the constraint ° { ^3 } <C 1, the integration over the Laplace transform is evaluated 
using the steepest descent method. It is reduced to 



^Low-Lying 

(W,V), 

1111 / 2 XV 14 

[6a {W 3 V} ) exp 



constant W V 



2Wy^V7 

11/ q \ 1/7 

where the saddle point of the stationary condition is found at 

lV7 



2/7 



i/r 



C 



(249) 



(250) 



The location of the saddle point must satisfy the energy constraint 
®{W 3 V} 1 (rat ~ \ 1/7 ^ 1 

6A^2 V 



6 A? 



1/4 



'6(2 j + l)iV c 



1/4 



7T Z 



[^3/4^1/4] <<; x 



' 6(2j + l)iV c n flavors \ 
J 

m < lOGeV (for 2 flavors). 



1/4 



6 TV 2 



«1, 



(251) 



In the realistic QCD, the quarks obey Fermi-Dirac statistics while the gluons obey Bose- 
Einstein statistics. The partition-like function for a gas of massless quarks and anti-quarks 
confined in a finite size cavity with the internal color structure SU(N C ) (or U(N C )) reads 

(2j + 1) 



a q q(R) 



7T Z 



(W 3 V)n^ V0T • tr / dxx 2 In 1 + R fun (<?)e x 



+ 



f 

Jo 



dxx In 



i + Kfun^K 



(252) 
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Similarly, the partition-like function for a gas of gluons is reduced to 

*<R) - - {23 + 1) 



(W 6 V) • tr f / dxx 2 In [1 - R adj (g)e 



(253) 



7T^ \J0 

It would be interested to note here that when a qq /( 3 < 0(l) cr ^j ca j and a g /( 3 < 0(l) cr ^j ca j, 
the color saddle points will be distributed uniformly over the entire circle. In this case the 
Maxwell-Boltzmann statistics can be produced as follows 



In 
- In 



i + RW^i-RfunG?)^ 

1-R adj (g)e- X ] » Radj^e-". 



(254) 



On the other hand, in the case a qq /( 3 > 0(l) cr j^j ca j and a g /( 3 > 0(1) critical the color 
saddle points become more dominant in a narrow domain around the origin. In this domain, 
the partition-like function is expanded around the color saddle points. The coefficients of 
the expansion around the saddle points for the quark and anti-quark partition-like function 
reads 

7tt 2 



a qq(R)\{e}=0 



d 2 a qq (R) 



09 



360 
1 



N c n RayoT (2j + l).(W 3 V), 



{6}=0 



6 



n fl 



avor 



(2j + l)-(W 3 V). 



(255) 



Moreover, the same thing can be calculated for the gluons 



a 9( R )l{0}=O 



7T 



d 2 a g (R) 



d{9 l - 9iY 



{9}=0 



-(N 2 -l)(2j + l)(W 3 V), 
= - l -(2j + l){W 3 V). 



(256) 



The coefficient of the quark, antiquark and gluon gas is reduced to 



a{w 3 v}(R0 = a qq (R) + a g (R). 



(257) 



The coefficient of the zero*' 1 approximation reads 



a {W [i V}(R) 



{9}=0 



= d qqg W 3 V, 



(258) 



where 



dqqg 



/7tt 2 tt 2 

W + (^"flavor + 5o(^-i: 



(259) 
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The density of states is approximated to the Gaussian-like function around the saddle points 
as follows 



P(ii)(W,V) = PHigh-LyingW V )> 

^3a {U /3y } (R)) /8 exp - (3a {H/ 3 y} (R)) 1 

2 qq ^ 2 2iV c Z ^ Z 1 V 



C qqg (W,V)exp 



9 s f 



where the pre-exponential coefficient reads 

1 1 



(3d qq - g W 3 v) 1/8 exp Uu mg W 3 V 



The quadratic terms in the exponential read 



i/4- 



A 



(2) 

qq 



' 1 \ 3/4 d\- q (R) 

sd q - qg w 3 v) de\ 



(3^ 5 ) 3/4 6 



«(2j + lHavor R 
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A (2) 
S 



3/4 



<9 2 a„(R) 



1 --(2j + i)(w 3 y 



(3<W 3/4 6 



. \ 1/4 



The color-singlet state for the density of states is projected as follows 



P(ii)(W,V) = C qqg (W,V) / dfi{g)exp 
= C qqg (W,V) [ C dfi(g)exp 



1 



N c 



L i=i 



2 2iV„ 



i=l j=l 



1 



^(A^ + 2iV c Af))E^ 



i=l 



The result of Eq. (12641) resembles Eq. (1167P and its evaluation leads to 

1 



P(n) 



(w,v) = Cq - qg ± (sd &g w 3 vf 



Ag ) + 2iV c A? ) 



exp 



(260) 



(261) 



(262) 



(263) 



(264) 



. \ V 4 ' 



(265) 



where the constant c g g g reads 



-'mi 



(2vr)^/ 2 n£ 1 fc! 



4AT c !(2tt)^- 1 /2 v ^FiV; 



(266) 
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In the terms of $^ and Eg. (1265 j) becomes 



4/„, Trr3Tr \V4" 

(267) 



where 



Af = (wV) 1/4 *f/(3d, 5! ,) 3 / 4 . (268) 

The Hagedorn density of states is given by P(n)(W, V) = PHigh-Lying^^' V)6(m— 4.BV) for 
the bags with the sharp surface boundary. It is found that Eq. ( 126 7ft is identical to the result 
given in Eq. ( 11 77ft . The constraint of the thermal running parameter for Eq. ( 12641) validity is 
given by 

* = W c (** + 2iYcA ? } ) - 5 ( Wh6re Ve-)critical = \) " ^ 269 ) 
Using the standard MIT bag model we get 



119 J 



mi^r^ + l) (" fla ™ r + 2Nc) m - 5- (270) 

where (W=V) 1/3 = 3 3/4 m./(4B 1 'i) an( j B 1 /! = 180 MeV. Eq.gg gives the minimum mass 
limit ^(27) m i n while the (pre-) critical mass threshold ^ cr jti ca l is determined by Eq. (I270p 
with the condition ^ cr jti ca l > m (//)min- The Hagedorn mass constraint is bounded from 
below by the (pre-) critical mass m > '^critical " Nonetheless, Eq. ( 12701) gives the following 
estimations 

m > "^critical = 2060MeV (massless flavors : ^fl avor = 1), 
m > ^ cr itical = 2315MeV (massless flavors : ^fl avor = 2), 
m > ^ cr itical = 2482MeV (massless flavors : ^fl avor = 3). (271) 

It is worth to note that Eq. ( 1225ft gives the constraint for the minimum mass limit for the 
Hagedorn threshold production where the (pre-)critical mass determined by Eq. (I270D must 
lie above this minimum mass limit ^ cr jti ca l > m (//)min- According to Eqs. (12261) and (12711) . 
the minimum mass threshold and the (pre-) critical mass for the two-flavor hadronic mass 
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spectrum are found m CT [^[ ca \ = 2315MeV > ^(/n m i n = 2029MeV. The experimental data 
book indicates that the maximum limit for the discrete mass spectrum is estimated to be 
around 2250-2300 MeV. The critical mass for the production of fireballs has to be just above 
the the known discrete mass spectrum. The present work predict that maximum limit for the 
two flavor discrete hadronic mass spectrum is 2315 MeV and above this limit the continuous 
Hagedorn spectrum is produced. 



X. THE GAS OF BAGS WITH VAN DER WAALS REPULSION 
A. Isobaric partition function 

The partition function for a gas of bags in the Maxwell-Boltzmann statistics reads 

1 N 

Z(T, V) = ET7yIl^(A), (272) 
n>o iV • i=l 

where 

Qi(A) = J d mi dv % J (v-Y,Vkj ^^p{k-m hVi )e- Ei l T Q (v - , (273) 

and Vi, rrii and Ei are the hadron's volume, mass and total energy, respectively. The fugacity 
A is the constraint of all possible charge conservations such as baryonic, strangeness, • • • etc 
and is given by 

A = { e «* . . .} = ■■■}. (274) 

The density of states given by Eq. fl272p consists all the discrete known mass spectrum 
particles and the continuous Hagedorn density of states as follows 

p(k]m i} Vi- ■ ■) = p {I) (m,v,T,A) + p (//) (m, v,T, A), (275) 

where the first term pm(m,v,T,A) corresponds the discrete low-lying hadronic mass spec- 
trum while the second term p(in(m,v,T,A) corresponds the continuous Hagedorn mass 
spectrum. The discrete low-lying density of states of the known mass spectrum particles 
reads 

Baryons 

p^(m,v,T,A) = D FD ( m '^ T ' Ai)5(m - m;)<% - v H ) 

i 

Mesons 

+ E D^(m,v,T,Ai)5(m - mi)5(v - v H ). (276) 

i 
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The sum runs over the baryon and meson mass spectra those are satisfying Fermi-Dirac and 
Bose-Einstein statistics, respectively. The terms Dt?d (m, v, T, A*) and -Dgg (m, v, T, Aj) 
are the degeneracies as well as the single-particle statistic ensemble functions for Fermi- and 
Bose- particles, respectively. The continuous part pnn(m,v,T,A) is the Hagedorn density 
of states and these states correspond the hadronic bubbles with relatively large hadronic 
masses could emerge as fireballs just above the highest mass of the known hadronic particle 
without strangeness m > 2.0GeV represented by the discrete low-lying mass spectrum. In 
the simplest approximation there is no reason to prefer the Fermi-Dirac or Bose-Einstein 
statistics for the exotic hadronic states such as Hagedorn states. The Hagedorn states are 
assumed to obey simply the classical Maxwell-Boltzmann statistics due to their relatively 
large masses. The grand canonical ensemble for a gas of non-interacting multi-particle 
species obeys the tensor product of their Fock spaces [381 ]. 

Z{T,V;A)= II Z ( T > V ) II Z (T, V) ]J Z(T, V) ■ 9 (v - $>V (277) 
Baryons Mesons Hagedorns ^ i ' 

In order to overcome the volume step function problem in Eqs. (l272p and (I273p . the isobaric 

ensemble trick is introduced. The isobaric partition function is calculated by taking the 



Laplace transformation of the grand partition function [35|, |36|, |37|, |38|, 139 ] 

Z(T,s;A) = / dVexp (—sV) Z(T, V; A), 
Jo 

= 1/ s - Jexp(-Vis)<p s (T;mi,Ai) , 

= V [ s -/Hadrons( T ' s )] » ( 278 ) 

where t>j is the hadron's Van der Waals excluded volume. It is reduced to a continuous Van 
der Waals variable v for the continuous mass spectrum. The isobaric function found in the 
denominator in the right hand side of Eq. (12781) is decomposed as follows 

/Hadrons( T ' s ) = / exp (—Vis) p(m, v, T, A)ip s {T; m*, A*), 

Ji 

= / ^exp(-t)is) p^(mi,v h T,A)(ps(T;mi,Ai) (discrete Low-Lying) 

i 

+ J dvdmexp (— vs) p(jj)(m, v, T, A)(p$(T; m, A) (continuous Hagedorn) 
= /Low-Lying ( T ' s ) + /Fireballs ( T ' s )- ( 279 ) 

The first term /Low-Lying s ) denotes the isobaric function consisting all the known 
hadronic mass spectrum particles and resonances as well as their antiparticles embedded 
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in the hot and dense medium. The masses of these particles are taken from the particle data 



book 



451 ] . The non-strange hadronic spectrum consists 76 mesons and 64 baryons. This 



isobaric function is given as follows 

Mesons 

4ow-lying( T ' s ) = S [^BE( T ; m " A ^) + ^BE( T ' m "^) 

i 

Baryons 



(280) 



The ensemble functions for single-particle species obeying the Maxwell-Boltzmann, Fermi- 
Dirac and Bose-Einstein statistics, respectively, read 



^ MB (T;mi,Ai) = (2^ + 1)^ J 



d 3 k 
W) 



g T 



^ FD (T;m J ,A J ) = (2J l + l) j ^ In [l + A.e^v 7 ^! 



and 



^ BE (T; m„ A,) = -(2J, + 1) / ^ In [l - A^V^I 



(281) 
(282) 

(283) 



where Jj is the quantum number stemming from the internal degrees of freedom (e.g. spin 
multiplicity). The effective chemical potentials read 



A* 



S A ): , A* = e~ vs Ai 



(284) 



where the isobaric pressure times Van der Waals volume (s ■ v) enters as an exponential 
pre-factor of the particle fugacity. The second term in Eq. (12791) corresponds the isobaric 
pressure of Hagedorn bubbles. The isobaric function for a gas of Hagedorn states reads 



•^Fireballs ( T ' s ) 



dv / dme vs p (II) (m,v, T, A)yj MB (T; m, A), 



(285) 



where the first integration is over the hadron's mass while the second one is over the hadron's 
excluded volume. The asymptotic behavior of the function V^MB^' m i A ) m ^ ne hmit of the 
large bag's mass m 3> T reads 

d 3 k 
~m 2 T~ 



/CI rC / \ 

^exp(-VA; 2 + mVT) 

K 2 (m/T), 



2tt 2 

mT\ 3 / 2 



/mi V 
\2n) 



-m/T 



(286) 
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By introducing Eq. fl286|) in Eq.([285]), Eq. fl285|) is reduced to 

/•co poo /rnT\ 3 ^ 

/Fireballs( T ' s ) = J dv J dme^'p^m, v, T, A) J e" m / T , (287) 

where p(n)(m, v, T, A) measures the mass spectral density and the volume fluctuation as 
we// of the continuous Hagedorn states. 

It is known from the property of the Laplace transformation that the asymptotic behavior 
of Z(T, V) as V — > oo is determined by the extreme right-hand singularity of Z(T, s; A) with 
respect to the isobaric variable s on the real axis. We denote this singularity point by s* . 
The Laplace parameter s* plays the role of the isobaric pressure. In the thermodynamic 
limit V — > oo, the pressure reads 

p = T Jim hnZ(T,V), 

V— >oo V 

= Ts*. (288) 

The isobaric partition function with the isobaric ensembles (T, s) is convenience for a system 
characterized by the external pressure p = Ts rather than the fixed volume V. The isobaric 
partition function has another singularity beside the first singular point So = s*. This 
singularity is determined by the nonlinear equation s* = /jjadrons s *) & yen by Eq. (12781) . 
The second pole s c is the singularity when the isobaric function /pireballs^' Sc ) diverges 
because the fireball's internal pressure exceeds the external pressure of the gas of hadrons 

/Fireballs ( T ' s c) -> oo. (289) 

The density of states for the Hagedorn states described by the standard MIT bag model 
with a sharp surface boundary reads 

p {II )(m,v) =5{m- ABv) Z (nm (W,v), (290) 

where W = m — Bv and the micro-canonical ensemble 



Z(II)qqg(W, v) = Cv~ N * /2 X \ 8 ) exp 



1 " x "4 .... 

-ax 3 ' f 



3 

(with (W, v)^(x = W/v, v)) , (291) 

is the micro-canonical ensemble for a quark and gluon gas projected on the color-singlet state 
and confined in a spherical cavity with a specific volume v. This micro-canonical ensemble 
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has been derived in sections IVIIII and IIXI Hence the isobaric pressure for the continuous 
spectrum of the hadronic bubbles becomes 

'mT\ 3 / 2 



roo roo / 777,7 \ ' 

/Fireballs( T > s ) = J vq dv dme- vs 6(m - ABv)Z {II)q - qg {W ) v) [—) e~ m ' T , 

= C'{T) ^ dvv~ N ' /2+3/2 exp [-v (s - s )} , (292) 



wo 



where the bag's isobaric internal pressure is given by 



s 



The integral pre-factor, which is independent on the volume fluctuation, reads 



(293) 



3iV:?+2 



/4 fi r r\ ^/^ 

C'{T) = C[3B]~\~s-J J • (294) 

In the case of QCD, we have iV c = 3, C*(T) = C'(T)\ N =3 and the isobaric pressure for 
Hagedorn's gas reads 

POO 

/Fireballs( T > s ) = C ^ T ) / dvv ~* ex P l~ v ( s ~ s o)] • (295) 

The pressure p = Ts* is typically calculated from the isobaric function extreme right singu- 
larity. The external isobaric pressure of Hadronic gas and the pressure of Hagedorn's gas, 
respectively, read 

s = ^Hadrons^' s )> 
^Fireballs = T ' /Fireballs( T ' s )' ( 296 ) 



where 



and 



/Fireballs^ «) = C*(T) (^V f° dzz-'e-, (297) 

Wo/ Jzn 



z =v (s- so) . (298) 
B. The order of phase transition 

In order to analyze the phase transition, it is useful to introduce the exponential integral 
function 

roo 

zl \ dzz~ n - x e~ z = z%r(-n,z ). (299) 
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This function diverges for n < as follows 

/•OO 

limzJ?/ tZ^- n_1 e-* = lim ^iY-n, z ) -> diverge. (300) 

^o^O z ->0 

In case n = it diverges logarithmically. On the contrary, it converges to 1/n for n > 0. 
Hence, the analysis to find the order of phase transition in the limit s — > s Q becomes a 
straightforward one. It is useful to adopt the following approximation, 

s = s Hadrons + C*{T) f°° d w - c -(-o ) 



1 

s^so 



s « C*(T) / dOT^e^' 8 ^"' 
Indeed, Eq.f l30ip leads to the following conclusions 



(301) 



23|, 



so 



for a > 1 : s — > finite (a possible phase transition), 

for a < 1 : s — > diverge (no phase transition). (302) 

In the case a < 1, the fireball pressure diverges in the limit s = so and subsequently the 
phase transition does not exist. On the other hand, in the case a > 1 the fireball isobaric 
pressure converges in the limit s = Sq, though it diverges as the bag's internal isobaric 
pressure exceeds the hadronic external isobaric pressure Sq > s and the phase transition to 
an explosive quark-gluon plasma takes place in the system. 

The order of the phase transition is determined by examining the continuation of the 
isobaric pressure outside and inside the Hagedorn's bag and its n-th derivative as well. Its 
derivative with respect to temperature (or any thermodynamical ensemble) reads 



s 



OO , POO 



' = (C*(T))' dvv- a e- v{s - so) -C*(T)(s' -s' ) / dvv- a+1 e~ v{s ~ so) . (303) 



VQ JV0 



The exponential integral in Eq. (13031) at the point of the phase transition leads to the following 
conditions 

finite value 



= 0, for 2 > a > 1 (a higher order phase transition), 

7^ 0, for a > 2 (a first order phase transition). (304) 

In the summary, first and second order phase transitions take place under the following 
constraints 

oo > a > 2, (first order phase transition), (305) 
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and 



2 > a > 1 + -, (second order phase transition), 



(306) 



respectively. The n-th order phase transition takes place whenever a takes the following 
value 



The Hagedorn bubbles considered in the present work are color-singlet states with an internal 
structure of a = 3. They may appear in RHIC and LHC as fireballs (i.e. the Hagedorn 
states). This means that the gas of fireballs (i.e. the Hagedorn phase) undergoes a first 
order phase transition to an explosive quark-gluon plasma. However, when the explosive 
quark-gluon plasma appears it expands rapidly. Nevertheless, in a context of an alternate 
scenario it is reasonable to assume that the colored bags could appear in the system before 
the appearance of an explosive quark-gluon plasma. The isobaric function can be written 
as follows 



/Fireballs( T . s )//colored-bags( T ? s )~/ dw^e'^ / dw^C^. (309) 



In this case the phase transition from the Hagedorn phase to the gas of colored bags is a 
first order phase transition. The surprise in this scenario is that the phase transition from 
the gas of colored bags to the explosive quark-gluon plasma is not possible, keeping in mind 
in this scenario we only consider a simple color group symmetry with no other associated 
symmetry. This leads to the conclusion that the explosive quark-gluon plasma must take 
place for the quark and gluon bags with some specific internal structure. 

On the other hand, it is also possible to think about a gas of colored bags with a spe- 
cific internal color-flavor correlation associated with additional configuration-space internal 
symmetry such as O(N) or SO(N) or even any other associated non-trivial effect. Then the 
generalization of the isobaric function for a gas of colored bags becomes 




(307) 



Finally, there is no phase transition for the bag of an internal structure of 



a < 1, (no phase transition). 



(308) 




(310) 
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with an internal structure | < a e < |. Then is this scenario, it is possible that the 
Hagedorn's gas f(T, ■s)Fireballs/^( T ' s ) exotic-bags undergoes a higher order phase transition 
to a gas of non-color-singlet bags. The resultant gas of non-color-singlet bags undergoes a 
higher order phase transition to an explosive quark-gluon plasma. 



XI. THE ROLE OF CHIRAL PHASE TRANSITION 



In this section the role of the chiral phase transition is considered in detail. A comprehen- 



sive review of the chiral fields can be found in Ref. [46]. In the context of the Gell-Mann Levy 
model, namely cx-model, the quark is assumed coupled to the chiral field of linear cx-model. 
The model can also be studied in the context of the Nambu Jona-Lasinio model in a similar 
way. The analysis of the chiral restoration in the context of the Nambu Jona-Lasinio model 
and color deconfinement will be considered elsewhere 471 ] . We can extend the Lagrangian 



density of the cr-model of Refs. 



48 



49] to include the flavor chemical potentials as follows 



C 



h^df, - m - gvMl^A^n - g a (a a X a + ry 5 A a 7r a ) + i- (8 B + 9 S • • •) 7 







+ \ [d^aa + d,7r a d^ a ] - l -F bllv Fr - U(a, vf), 



(311) 



where the axial scalar and axial pseudo-scalar field potential is given by 



^■6 I 



+ 



12 



K 



1 



tr F [(a a A a + in a \ a )(cr a \ a - irr a \ a )] 2 - - [tr F (a a A a + i7r a A a )(cr a A a - i7T a \ a )] 2 



— [det (a a X a + i7r a A a ) + det (a a X a + m a K)] ■ 
The non-Abelian color-gluon field reads 



(312) 



Fbpuv — d^Abu — dyA^ — gy m fovb" Aw nAwv. 



(313) 



The quark satisfies the color fundamental representation defined by the following group 
transformation ^f un = exp {i9 a T a ) where r a is a set of the fundamental symmetric group 
SU (N c ) generators. On the other hand, the gluon satisfies the adjoint group representation 
defined by the transformation = exp (i(fi a T a ) where T a is the set of adjoint generators of 
the same group that generates the fundamental generators. Furthermore, the adjoint color 
parameters 4> a are related to the fundamental ones by the relation <p a oc (0j — 9j) where the 
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index a runs over 1, • • • , (N 2 — 1) while i,j indices run over 1, • • • , N c . The adjoint and 
fundamental indices are related by a = (ij). The set of generators that commute with the 
Hamiltonian is retained with the correspondence conservative parameters set 0j. The details 
will be given elsewhere \^\ - The last term in the Lagrangian 

Vui = — — [det(a a \ a + m a K) + det(a a A a - iTT a X a )} , (314) 

is t^ ax j a i (1) symmetry breaking term. For two flavors Nf = 2, we have r/^ = 0. The model 
can be simplified furthermore by retaining only one scalar a field and three pseudo-scalar 
iii fields and the constant k 2 = of the axial symmetry breaking term. The easiest way to 
analyze the chiral phase is to adopt the mean-field approximation by replacing the a and 7Tj 
fields by the their expectation values or condensations a ps (a) and m ~ (71*), respectively. 
The trivial example is the set for the symmetric nuclear matter (7r») =0. Finally, we neglect 
the interaction between the gluons and quarks by setting the non-Abelian coupling constant 



to qym = 0. The small coupling qym can be treated perturbatively [47]. 
The effective chiral Lagrangian reads, 

C = C q - q + C g + C a . (315) 

The quark's term becomes 

Cqg « q 



q 



q, (316) 



where chiral constituent quark mass is given by m*((a)) = mo+ g a (&), while the conservative 
charge parameter such as the baryonic charge is defined by iQb — ^fr- Hereinafter, we shall 
neglect the chemical potential /i^ as far we are interested in the phase transition along 
temperature axis for the diluted and hot nuclear matter. The hot hadronic matter at zero 
baryonic density simplifies the calculations drastically. The gluon-Lagrangian part in the 
limit of zero coupling constant, gyu = 0, is reduced to 

C g « -\f%F^\ (317) 

where the color-gluon field becomes 

F; v = dpAZ - QyAl. (318) 
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There is an additional term in the Lagrangian due the chiral interaction. This chiral- 
Lagrangian term is responsible to the chiral restoration phase transition at high temperature. 
The chiral-Lagrangian reads 

1 



-<r ~ jd^a- U(cr), 



-u ({*))■ 



(319) 



Hereinafter, the mean field brackets are discarded and redefined as (a) = a. 

The partition function for a gas of quarks and gluons and chiral source is given by a 
tensor product of the quark, the gluon and the chiral field Fock spaces 



Z qqga (P,V;a) = Z q - q (P,V;a)-Z g (P,V)-Z a (P,V;a). 



(320) 



However, the chiral potential is to be subtracted from the quark and gluon bag energy 
when the micro-canonical ensemble is computed. The whole idea is that the hadronic bags 
are embedded in the chiral field background and the constituent quarks are coupled to the 
external chiral field. 

The color-singlet state is projected as follows 

Z SingleW V ~> a ) = J dfi(g)Z gqga ([3,V; a), 

' dn(g)Z qq (P, V; a) ■ Z g ((3, V) ■ Z a ((3, V; a), 



J dLi(g)Z qq ((3,V;a)-Z g (f3,V) 



Z a (P,V;a). 



(321) 



The quark-antiquark partition function with an internal color degree of freedom is given by 

Vd 3 p 



Z gq ((3,V;a) = exp 
•In 



i + Kun^ + Rfun^K 



{2nf 



(322) 



The gluon partition function is not modified by the chiral field and it reads 



Z g ((3,V) = exp 



(2 J + l) f Vff ln 



l-R adj (s)e"V^ 



(323) 



7V C 2 - 1 " C J (2tt) 3 

where the gluon mass remains m g = in the hadronic phase. Finally, the a-energy partition 
function is calculated from the a-potential as follows 



(324) 
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The micro-canonical ensemble is calculated by finding the inverse Laplace transform of the 
mixed-canonical ensemble 

1 r/3 c +oo 

Z{W > V]0) = 2^L d ^ WZ SmgletA^ V ^)- ( 325 ) 

Z7T6 J p c — 'too ° 

The a-potential background is to be subtracted from the bag energy in order to scale the 
quark and gluon bag's energy correctly. The cr-field is the scalar chiral field interaction 
among the bags and it is coupled to the constituent quarks. The quark and gluon bag's 
energy W is scaled and re-defined as follows 

W = W -VU a (a). (326) 

Hence, the micro-canonical ensemble of the chiral quark and gluon bag is calculated as 
follows, 

Z{W > V;a) = 2^n jLoo d ^ W+ " VU ^ Z Sin g \eU (ft V; *) > 

= ^/ ft - to> d ^ ro - WZ Singlew(A^-)- (327) 

The resultant micro-canonical ensemble given by Eq. (1327p is equivalent to start from the 
beginning with a mixed-canonical ensemble as a tensor product of the Fock spaces of chiral 
quark and gluon projected in the color-singlet state as follows 

Z(W, V; a) = — / df3e^ w Z Sia let ((3, V; a) , (328) 

where W is the quark and gluon bag energy while the chiral color-singlet mixed-canonical 
partition function is given by 

^Singlet (Z 3 ' V 'i a ) = J M9)Z q q((3, V; a) ■ Z g ((3), 

= Jdfi(g)Z q - qg ((3,V;a). (329) 

It is demonstrated in sections IVHI and IVIIII that the non-chiral hadronic phase under- 
goes a third order phase transition from a hadronic gas dominated by the low-lying mass 
spectrum to another hadronic gas dominated by the continuous Hagedorn mass spectrum. 
It is interpreted that low-lying mass spectrum corresponds to the discrete hadronic mass 
spectrum found experimentally and consisting meson, baryons and any exotic hadronic mass 
spectrum could be found in the data book[45j. On the other hand, the Hagedorn states are 
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the highly excited metastable hadronic states those are produced just above the discrete low- 
lying mass spectrum. The density of states for the continuous Hagedorn mass spectrum is 
calculated from the micro-canonical ensemble. In the calculation of the color-singlet canon- 
ical ensemble, the multi-integrations over the colors are preformed. The standard procedure 
is to adopt the saddle points approximation. In the limit of chiral Hagedorn states, the 
quark and gluon canonical ensemble is approximated by expanding the exponential around 
the color saddle points up to the quadratic term similar to the same procedure preformed 
to the non-chiral Hagedorn states. The quadratic expansion of the chiral quark ensemble 
given by Eq. ( 13221) with respect to the color saddle points reads, 

In Z q - q ((3, V, a) = $W - \*% £ 0* (330) 



i=i 

The coefficients of the zeroth and the quadratic terms read, respectively, 



= 2N I ' - - 2 In 



Vd 3 p 



V2N r . 



and 



(3 3 IT 2 

frnfaj) , (331) 



$ ( 2 ) = m l_ [Vd 3 P 2e-V^^ 



y 2 f°° 2 e ~\/^ 2 +/3 2 m% 2 (<r) 



(3 3 7T 2 Jo 



f 1 _|_ e -^+(Pm% 2 {a) 



2 ! 



— $ ( - } ((3m* (a) 



V 



^%((3m* q {a)), (332) 



where m* (a) = (mo + (ja(y\ The first order term is not needed in the expansion due to the 
saddle point constraint of extremization. 

The same thing can be performed for the gluon-part given by Eq. fl323p . The only differ- 
ence is that, in the gluon ensemble, the expansion is carried over the adjoint color variables 
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rather the fundamental ones in the chiral quark- ant iquark ensemble. The gluonic quadratic 
color expansion is identical to that found in the non-chiral Hagedorn as follows 

i N c N c 

ln^ 9 (/3) = $f - r*f E E (*» ~ 6 m ) 2 , (333) 



2 n=l m=l 



where 



$f = where $f = 2(JV? - 1)^, 

$f = |^ 2) , where $f = ±. (334) 

The resultant mixed-canonical ensemble that generates the density of states for the color- 
singlet chiral Hagedorn Eq. (I329[) is approximated to 

Z Singlet(^'^ a ) ~ J d K9) S Q,ddle Z i9(^ V ^ a ) Z 9(^ V )^ ( 335 ) 

where / d[i{g) SSu ddle * s ^ ne hivariance measure given by the Gaussian-like saddle point 
method and Eq. (1221) . After evaluating the multi- integrations over the color variables us- 
ing the saddle points approximation, the chiral color-singlet mixed-canonical ensemble is 
reduced to 



*<£ ' (/9m£)+*< 0) 



.T— J. ,336, 

where the pre-factor normalization coefficient M q q g is identical to the non-chiral mixed- 
canonical ensemble given by Eq. (1169l) . 

The chiral color-singlet micro-canonical ensemble is found by calculating the inverse 
Laplace transform of the mixed-canonical ensemble as follows, 



Z (II)q - qg (W,V;a) = dp ( ~ (2) 1/V (2U - 



x exp 



qq, 

— (&°l(3m*)+& 0) 



(337) 



The inverse Laplace transform is evaluated using the steepest descent method. The Laplace 
stationary- saddle point is determined by extremizing the exponential under the integral in 
Eq. (133 7p with respect to the Laplace transform variable (3. Unfortunately, the solution of the 
Laplace-saddle point is found a transcendental one and can not be written in a closed form. 
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However, we are concerned by the solution near the chiral restoration phase transition where 
the constituent quark masses approach their current ones. In the limit of small constituent 
quark masses such as light flavors, the Laplace saddle point is found by iteration. In this 
kind of transcendental problem, the solution converges rapidly. The saddle point solution is 
found as follows, 



(i) 



ft 



% (*Sl(/K « °) + *W 



A) » ft 



(In) 



and 



(338) 



where I n is the number of iteration. In the present kind of transcendental equation, we 
assume that the value (3q = (3^ (i.e. I n = 2 and f3 = (3^) is a sufficient approximation in the 
limit below but close to the chiral restoration phase transition with small constituent quark 
masses like up and down flavors. This approximation is also justified for the strangeness. 
The micro-canonical ensemble takes the following form, 



Z (n)q - q9 (W,V;a) 



1 



2V27T 



I (3iVg+2) 

" W 8 exp 



1 w*'*vw 



(339) 



for the chiral color-singlet bag of quarks and gluons. The pre-exponential coefficient N$ is 
dimensionless and reads 



(340) 



(2iV c $f + $gt(/3m*))" 



The neutralness of the pre-factor coefficient N$ does not generate a power function with 
respect to the bag's energy or volume. In this sense, its variation with respect to the 
chiral constituent quark mass does not play a significant role in the chiral restoration phase 
transition and subsequently this variation is neglected. Therefore, it is reasonable to ignore 
the chiral effect in N$ and to approximate it to the following: 



(*SU(o)) : 



(2iV c $ 



f + a. 



(2) 
qqa 



(0))' 
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(2 • 2N C • 7tt 2 /720 + 2 • (iV c 2 - l)7r 2 /90) (3iVc2 ~ 2)/8 



(341) 



(2iV c /6 + l/3) { ^ 1)/2 

In the contrary, the exponential plays a significant role in the chiral phase transition. It 
is decisive to determine the scalar chiral a-mean field and the order of the chiral phase 
transition. It must be written explicitly as a function of the chiral constituent quark mass 
as follows 



(342) 



The explicit expression with respect to the chiral constituent quark mass reads 



qqa 



{(5ml) 



7T" 

3tt 2 



dxx In 



dx- 



x 



1 



1 + e 



(343) 



^x 2 + ((3m* a ) 2 

The variation of the exponential with respect to the scalar chiral cr-mean field is essential 
for calculating the scalar chiral cr-field. This can be done by extremizing the isobaric pres- 
sure. The explicit variation of $^(/3m*) with respect to the scalar chiral cr-field takes the 



following form, 
5 



5a 



2N C r 

7T 2 JO 



dxx 



f3 2 m* 



V 



;V /x2+(/3mJ) 2 + l 



x 2 + (/3m*) 



J 



5a 



(344) 



The variable transformations x = W/V and v = V transform the independent variable set as 
follows {W, V} — > {x,v}. However, this transformation set simplifies the analysis of phase 
transition in terms of the bag's energy density and volume rather than the bag's energy and 
volume. Nonetheless, according to the standard MIT bag model, the bag's mass is given by 
subtracting the bag's energy constant as follows W = m — Bv, where BV is interpreted as 
the bag's vacuum energy. The micro-canonical ensemble as a function of two independent 
variables {x, v} is displayed as follows 

"4 



<(n)qqg{X,V,a) 



1 . nI 
Nqq9 N S V 2 X 



cxp 



{^qi^oy^v 



(345) 



2V27T 

The Hagedorn's density of states in the context of MIT bag model with a sharp surface 
boundary becomes, 



PBags = S (m - 4Bv) Z {II)q g g (x,v;a), 
= ^8(x- 3B) Z( II)q q g (x, v; a). 



(346) 
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The vantage of point adopting the standard bag model with a sharp surface boundary 
condition is simplifying the model drastically. However, in the realistic physical situations, 
it is expected that the bag surface boundary has an extended surface in the extreme hot bath 
rather than a sharp one. Furthermore, it would be reasonable to think that the bag's surface 
boundary is distorted and becomes a fuzzy one. It is argued that bag's surface boundary 
is important in the determination of the order of phase transition and the existence of 
the tri-critical point in the phase transition diagram 0, [|3(3]. Nonetheless, the surface 
distortion may cause a mechanical instability. Such a mechanical instability is ignored here. 
The Hagedorn's isobaric pressure is calculated by integrating the degenerate single-particle 
ensemble function over the mass and volume distribution function. It is reasonable to 
assume that the gas of Hagedorn states is satisfying the Maxwell-Boltzmann statistics in 
order to simplify the model analysis. The Hagedorn isobaric function is integrated over the 
continuous mass and volume variables as follows 

/*oo /*oo 

/Fireballs ( T ' s ) = / / dvdme~ vs p B ^{m, v)(p MB (T; m). (347) 

As far the initial Hagedorn mass is relatively heavy, it is adequate to adopt the following 
approximation, 

/oo poo 
J mo dvdme ~ VS PBags( m i v ){-^) e > 

= r r dvdxe-™8(x - 3B)Z (II)qqg (x, v; a) f ( x + B ) vT \ ^ e -{*+B)v/T 

Jv() Jx() \ Z7T / 

(348) 

The integration over the ^-function eliminates the energy density integration and reduces 

the isobaric function to an integration only over the Hagedorn's excluded volume as follows 

/■oo /4 R7^\ 3 / 2 

/FireballsCM = / q d™"' < =3B (^) v^e~W T >. (349) 

In order to simplify the analysis, the chiral isobaric function for the Hagedorn gas is written 
in short as follows 

/FireballsCM = C '( T ) H dw-^e^-^, (350) 

Jvn 



where the integral's pre-factor function reads 



1 3ivg+2 /4BT\ 3/2 

C\T) = ^j=M q - qg N~{Wy^ (— j . (351) 
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It is worth to note that Eg. (13501) differs from Eg. (1295[) where the former equation includes 
the scalar chiral cr-field unlike the later one. Furthermore, the Hagedorn bag's internal 
isobaric pressure is given by 

^Internal /-^ = s °' 

= \ {*%*<Mt)) k (35) 3/4 - 4B/T. (352) 

On the other hand, the chiral low-laying mass spectrum is realized as the discrete hadronic 



mass spectrum found experimentally and displayed in the data book [45j . The scalar chiral 
cr-mean field potential interacts and couples to the hadrons and generates the hadron masses. 
At highly extreme conditions, the chiral symmetry will be restored and the hadrons with light 
flavors will dissolve to massless states while the scalar chiral cx-potential vanishes. Under 
the above assumption the density of states for the low-lying mass spectrum is generalized 
to take into consideration the chiral discrete mass spectrum as follows, 

Baryons 

PLow-Lying( m > u > T > A ) = £ ( 2 J B + X ) D FI)( A b) S ( m ~ m B^) 6 ( v ~ v h) 

B 

Mesons 

+ E (^U + 1 )^E( A M) s ( m - m M^) 6 ^- VH ')- 



M 



(353) 



The sums run over the baryon and the meson mass spectra those are satisfying Fermi-Dirac 
and Bose-Einstein statistics, respectively. The coupling of the scalar chiral cr-field to the 
constituent quark generates the hadron's mass as follows 

M Hadron = 9H«o + n q m q0) 

« n q ■ 4^ ut • a, (354) 

where n q = 2, 3 for meson and baryon, respectively, and m q o is the quark current mass and 
it is usually massless for light flavors m q o ~ 0. In general, the coupling constants gua or 
Qqcr^ are determined by the phenomenology in order to fit the experimental hadronic mass 
spectrum. However, the above procedure is boring and produces a tremendous number of 
fitting parameters. In order to reduce the number of phenomenological parameters gn a and 
to include the effect of the discrete hadronic mass spectrum, we assume the effective discrete 
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hadronic mass spectra for baryons and mesons are, respectively, generated by 

m* B (a) = (M B - M Nudeon ) + 3g qa a, 

m* M (a) = (m m - -M Nucleon j + 2g qa cr. (355) 

In the computational calculations, a special consideration is given for light mesons such as 
pions as far they are replaced by the scalar a-mean field. Furthermore, the u meson is 
replaced by the vector mean field while the other non-strange hadrons are left as hadron 
particles. The generalization of the isobaric function of the discrete low-lying hadronic mass 
spectrum, namely the chiral low-lying isobaric function, is given by 

Baryons 

/Low-Lying ( T ' s ) = 12 [v?FD ( T > s > m B( a )' A b) + ^FD ( T > s > m h( a )' ^b)] 

B 

Mesons 

+ 12 [^BE ( T ' s ; m M( a )'%) +^BE ( T ' s ' m M( a )'^M) • 
M 

(356) 

The ensemble functions for single-particle species obeying the Fermi-Dirac and Bose-Einstein 
statistics read, respectively, 



V?FD {T, s; m* B (a), Ag) = (2 J B + l)J^\n 



1 + A B e~ vs e 



, (357) 



and 



V?BE (T, s; m* M {o), A* M )) = - (2 J M + l) J -0^ In 



_ 1 / fc 2 +m .2 ( CT ) 

l-A M e-^e T V M y 



.(358) 



The pressure p = Ts* is determined by finding the isobaric extreme right-hand singularity 
as follows, 

s * = /Fireballs ( T ' s *) + /Low-Lying ( T > s *) + /-( T > s *)> 

= /Fireballs( T ' s *) + 4ow-Lying( T ' s *) ~ U W)' ( 359 ) 

where the isobaric term of the scalar chiral a-field is given by f a (T, s*) = —U(a). 

The value of the scalar chiral a-mean field is determined by extremizing the chiral isobaric 
pressure as follows 

^-s* = 0. (360) 
da 
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The chiral variation of the Hagedorn's isobaric function reads 

^/Fireballs^ O = C{T) £° d W -^ e -(-o)^ + (•••) ^. (361) 
The chiral variation of the bag's internal pressure reads, 

^Internal = T ^Q_ / 362 n 
da da ' 

where 

^ = l {3B 0^ K) f'\»M<l. (363) 

The asymptotic approximation of the chiral variation of the Hagedorn's isobaric function 
reads 



^/Fireballs( T ' O « ~K / q ^"^e^"-) + (• ■ •) ^- (364) 



The variation of the chiral isobaric function of the discrete low-lying mass spectrum with 
respect to the scalar chiral cr-field reads 

r, Baryons 

75Z WLyingC^ s *) = E [4d ( T ' s " m B^)> A b) + 41) ( T ' s " ™B^)> ^B 
r B 

Mesons 

+ £ [v^BE ( T ' s * ;m M( a )' A M) +^BE( T ' s;m M( a )'^M 
M 

+ (-)^. (365) 
The explicit expressions used in the above equations are displayed as follows 

^FD( T ^*^B( ff )' A B) = 4D( T ' s " m B^)' A B) + (---)^, 

V%(T,s*;m h (a),A h ) = - ( 2 J B + i)g v f™L. V ^ 66) 



AgV s e V B + 1 



and 



43,67) 



2vr 2 



A -i e , s . e V fc2+m M (CT) 
M 
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where M a (a) is interpreted as the effective mass of the scalar chiral cr-field particle. The 
dots term is the remaining terms those do not affect our final interpretation and conclusion 
in the chiral restoration phase transition. 

The scalar chiral a-mean field condensate is determined by extremizing the total isobaric 
pressure (i.e. grand potential) in the following way: 



It has been pointed out that the hadronic phase exhibits the Gross- Witten point in a way 
that the gas of the discrete low-lying mass spectrum undergoes a third order phase transition 
from the discrete low-lying hadronic mass spectrum phase to a hadronic phase dominated 
by the continuous Hagedorn mass spectrum. This leaves three possibilities in order to 
study the chiral phase transition. In the case (I), we have a mixed low-lying and Hagedorn 
phases near the Gross- Witten point. In this case, the Gross- Witten and chiral restoration 
phase transitions overlap around the Gross- Witten point, although the chiral restoration 
phase transition is a cross-over one unlike the third order Gross- Witten phase transition. 
Moreover, the both transitions take place below and far away from the deconfinement phase 
transition to an explosive quark-gluon plasma. In the second case (II), the chiral restoration 
phase transition takes place in the hadronic phase dominated by the discrete low-lying 
hadronic mass spectrum particles. In this scenario (i.e. case (II)) the chiral restoration 
phase transition is a smooth cross-over and is located below the Gross- Witten point and far 
away from the phase transition to an explosive quark-gluon plasma. In this case (II), the 
scalar chiral cr-mean field is determined by solving the following equation of extremization, 




(369) 



The Eq. (13691) is approximated and simplified to 




(370) 




— r 

da - 




(371) 
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This equation produces a smooth cross-over chiral restoration phase transition where the 
scalar chiral a-mean field decreases smoothly. This class of equation has been studied 
extensively in context of the Walecka model where the Van der Waals effect is neglected. 
When the bag's excluded volume is ignored, we get the standard a-model as follows 

d 



^Low-Lying (-^> s 



2(/B + 1 )^£/ 
B 



dkk' 



m*(cr) 



M 



2tt 2 



dkk : 



fc 2 +mjg(cr) 



m*(a) 



2vr 2 



(372) 



where the fugacity are set to Ag = 1, Avj = 1 as far the present analysis focuses in the region 
near zero baryonic chemical potential along the temperature axis. In this case, the chiral 
restoration phase transition takes place before the Gross- Witten point in contrary to the 
case (I) where in the case (I) the chiral restoration phase transition overlaps the the Gross- 
Witten point and is located far away from the deconfinement phase transition to an explosive 
quark-gluon plasma. Hence, the case (II) sounds to fail to produce the continuous Hagedron 
threshold since it generates massless hadron masses before the continuous Hagedorn mass 
spectrum is reached. Since in this scenario, the Hagedorn mass threshold is much higher than 
the maximum discrete mass spectrum, the continuous Hagedorn states production is denied. 
The self-consistent phase transition from the hadronic phase to an quark-gluon plasma is 
not possible in the context of case (II) because the gas of bags will never be produced unlike 
the cases (I) and (III). Hence, the Gibbs construction becomes essential in case (II) if we 
want to study the phase transition to quark-gluon plasma. It is hard to image an explosive 
quark-gluon plasma phase transition in this case. However, the case (I) sounds to be more 
physical than the case (II) since in the former case the chiral restoration phase transition 
overlaps the third order Gross- Witten point. This overlapping mechanism virtually make 
the Hagedorn phase difficult to be detected directly experimentally 51]. Furthermore, in 
this case the chiral restoration phase transition is supposed to take place prior to the color- 
deconfinement phase transition to an explosive quark-gluon plasma. Technically speaking, 
it will be hard to distinguish the color-deconfinement from the chiral restoration phase 
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transition. Nonetheless, the indirect detection and ana 



ysis can distinguish the explosive 
quark-gluon plasma from the static hadronic phase [52|, [53] . Finally, the third case (III) is the 
most important case. Indeed, in this case the chiral phase transition takes place in a hadronic 
phase dominated by the Hagedorn's gas. The chiral restoration phase transition takes place 
above the Gross- Witten point and utmost simultaneous the point of the deconfinement phase 
transition to an explosive quark-gluon plasma. The scalar chiral cx-mean field is determined 
by extremizing the isobaric pressure as follows 

\-L « Ml(a)a, 

a 

8 r i 

= -fo [/Fireballs ( T > S *)J ■ ( 373 ) 

Hence, the hadronic matter dominated by the fireballs (i.e. the Hagedorn states) leads to 
the following equation 

Ml(a)o = (•••)<•/ dvv- a+l e-' u{s - So) . (374) 

JVQ 

Eq. (I374p has very important characteristics which has been discussed in details in section [X] 



and Ref. [28(. At the point just below the phase transition from the Hagedorn phase to an 
explosive quark-gluon plasma, it is shown that Eq. (13741) has the following characteristic 
properties, 



MAa)a oc m* • lim (s - s ) a T (-a + 2,v (s - s )) , 

y s— »so 

~ finite for a > 2, 

~ oo for a < 2 -> a = 0, (375) 

where T(—n,x) is the exponential integral function of order n where n > 0. The integral 
on the right hand side of Eq. f)374p diverges whenever the bag's internal isobaric pressure Sq 
exceeds the external one s. The phase transition from the hadronic phase to an explosive 
quark-gluon plasma takes place in the limit (s — sq) — > when the bag's internal isobaric 
pressure so reaches the external isobaric hadronic pressure s from below. The solution 
a = is the only non-trivial solution for the scalar chiral cx-mean field condensate when the 
explosive quark-gluon plasma takes place subsequent the point (s — s ) = CT. The abrupt 
vanishing of the scalar chiral cx-mean field indicates evidently that the chiral symmetry 
restoration phase transition takes place below or utmost at the same point of the phase 



transition to an explosive quark-gluon plasma. Actually, the point s = so decides the 
order of chiral phase transition where the chiral phase transition takes place simultaneously 
with the deconfinement phase transition to an explosive quark-gluon plasma. However, 
if the cross-over chiral restoration phase transition takes place prior to the deconfinement 
phase transition then it takes place far away from the explosive quark-gluon plasma phase 
transition. 

The same analysis done in section |X] to determining the order of phase transition can be 
carried out for the chiral symmetry restoration phase transition. As the isobaric ensemble 
reaches the limit (s — so) = + just below the point of the phase transition, the scalar chiral 
u-mean field solution is finite for the gas of bags with an internal structure of a > 2 while 
it must vanish for a < 2 in order to get the physical solution. This means that for the 
gas of bags with internal structures of a < 2 the chiral restoration phase transition takes 
place below the point of the deconfinement phase transition to an explosive quark-gluon 
plasma. In this case, the order of the phase transition will be a rapid cross-over transition. 
Furthermore, the chiral restoration phase transition must be completed far away the point of 
the phase transition to an explosive quark-gluon plasma is reached. If the smooth cross-over 
phase transition is not completed prior the deconfinement point, then the chiral restoration 
is prohibited for a < 2. On the other hand, the case is completely different for the gas of 
bags with the internal structure of a > 2. In this case, the gas of bags undergoes the chiral 
restoration phase transition simultaneously at the same point of the explosive quark-gluon 
plasma phase transition. The order of the chiral restoration phase transition is found of a 
second order for a = 3 while of a higher order for 3 > a > 2. Finally, the chiral phase 
transition is of the first order for the gas of bags with an internal structure of a > 3. 

In the context of a model consisting hadronic matter dominated by the fireballs with an 
internal structure of a = 3 (i.e. the color-singlet bag), the scenario of the phase transition is 
summarized as follows: The gas of the discrete low-lying hadronic mass spectrum particles 
undergoes a third order phase transition to the Hagedorn phase at the Gross- Witten point. 
The chiral restoration phase transition likely takes place in the hadronic matter that is dom- 
inated by the continuous Hagedorn states. The chiral restoration phase transition is located 
above the critical Gross- Witten point and it takes place simultaneously with the point of the 
phase transition to an explosive quark-gluon plasma. The chiral restoration phase transition 
is of the second order and it takes place just below or utmost simultaneously the point of 
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the phase transition to an explosive quark-gluon plasma. Finally, the Hagedorn phase un- 
dergoes a first order phase transition to an explosive quark-gluon plasma. Nonetheless, this 
scenario does not exclude the possibility for the sharp cross-over chiral restoration phase 
transition to take place subsequent Gross- Witten point but far away from the point of the 
phase transition to an explosive quark-gluon plasma. However, in the case that the smooth 
cross-over chiral restoration phase transition is not completed yet in the Hagedorn phase, 
then the second order chiral restoration phase transition must take place simultaneously 
with the point of phase transition to an explosive quark-gluon plasma. Indeed, in this sce- 
nario it seems likely that the chiral restoration phase transition takes place in the continuous 
Hagedorn spectrum phase rather than in the discrete low-lying mass spectrum phase. The 
continuity of the Hagedorn mass spectrum makes the Hagedorn phase domain narrow in the 
hot bath. The narrowness of the Hagedorn phase range besides being the restoration phase 
transition lain above Gross- Witten point and taken place utmost simultaneously with the 
deconfinement phase transition makes the continuous Hagedorn states difficult to be dis- 
tinguished from the quark-gluon plasma. It makes also more hard to distinguish the chiral 
restoration phase transition from the deconfmement phase transition. For example it has 



been argued that the heavy Hagedorn states have large widths [5l|- Nevertheless, this sce- 
nario, however, is inconclusive although it explains the existence of the Gross- Witten point, 
the chiral restoration, and deconfmement and even the existence of the tri-critical point for 
some class of bag's internal structures. The existence of the tri-critical point is essential in 
the multi-processes scenario. 

The other alternative scenarios for a hadronic gas dominated by the fireballs can also 
be imagined. The fireball's internal structure a varies under extreme conditions due to the 
color-flavor correlations and the modification of the bag's volume fluctuations. In the case 
that the bag's internal structure is of order a > 3, then both the chiral restoration and the 
color deconfinement are of first order phase transitions if they are taken place simultaneously. 
In contrary, the chiral restoration phase transition is found of a higher order for the gas of 
bags with the internal structure of 3 > a > 2. The phase transition to an explosive quark- 
gluon plasma is of the first order for a > 2 while of a higher order for 2 > a > 1. 

In the class of scenarios where the chiral restoration and deconfinement phase transitions 
overlap, the n-th order chiral restoration phase transition takes place simultaneously with 
the first order phase transition to an explosive quark-gluon plasma. Nonetheless, the Gross- 
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Witten point remains a third order phase transition point. However, the simultaneous 
phase transition to the chiral restoration and to an explosive quark-gluon plasma is not 
conclusive. If the fireball has the internal structure of 2 > a > 1, then the deconfinement 
phase transition is of n-th order while the chiral phase transition must take place far away 
from the point of phase transition to an explosive quark-gluon plasma. Moreover, the gas 
of hadronic bags with the internal structure of a < 1 never undergo the phase transition to 
an explosive quark-gluon plasma. It is interesting to note that much of the above scenarios 
are precursors that the smooth cross-over chiral restoration phase transition is taken place 
below and far away from the point of the phase transition to an explosive quark-gluon plasma 
and moreover as the chiral restoration phase transition approaches the point of the phase 
transition to quark-gluon plasma the cross-over restoration phase transition becomes more 
sharp and rapid. 



XII. CONCLUSIONS 



In order to study the chiral restoration and deconfinement phase transitions in QCD, the 
hadronic density of states must be known for the entire energy domain below the point of the 
phase transition to the true deconfined quark-gluon plasma. The theoretical procedure to 
find the hadronic density of states is carried out by computing the micro-canonical ensemble. 
The canonical ensemble for the color-singlet bag of constituent particles with the underlying 
symmetric group SU(N C ) (and U(N C )) has been considered in detail. The color structure 
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34j or even the extended gauge field theories such 



existence of the tri-critical Gross- Witten point. This point is thought to play a significant 
role in the phase transition mechanism. It has been derived using the spectral density of 
color eigenvalues method originally introduced by Brezin et. al. [8|. Recently, the Gross- 
Witten critical point and the deconfinement mechanism has been reviewed in the context of 
the Brezin et. al. js] method to study the quark-gluon plasma in QCD 0, 1(3, 11, 13, Q, 3] 



and the black hole formation in AdS/CFT 22 



24. 25. 



26|. 



In the present work, we have introduced an alternate and simpler method to derive the 
color-singlet canonical ensemble for the asymptotic large thermal running parameter X/N^. 
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Furthermore, we have demonstrated in details how to locate the Gross- Witten critical point 
for the phase transition. This novel method suits the realistic and complicate physical sit- 
uations. For the small thermal running parameter \/N%, the saddle points are distributed 
uniformly over the entire color circle range \9i\ < n. Since the saddle points are distributed 
uniformly over the entire range, the Vandermonde determinant contributes to the action 
as an additional effective potential term. The integral of the resultant ensemble has been 
evaluated trivially. However, this procedure fails when the saddle points congregate around 
the origin rather than distribute uniformly over the entire color circle range \9i\ < it and 
the Vandermonde effective potential develops a virtual singularity. In this case, a further 
consideration must be taken into account in order to regulate the action involving the Van- 
dermonde effective potential and finally to evaluate the canonical ensemble correctly. There- 
fore, this behavior indicates that the solution changes its analytical function characteristic 
and a subsequent phase transition takes place. The solution of the asymptotic large A/iV? 
is evaluated using the Gaussian-like saddle points method. The action is expanded around 
the stationary Fourier color variables where only the quadratic terms are retained while the 
Vandermonde determinant is regulated in a nontrivial way. Fortunately, in this procedure 
the stationary Fourier color points are found dominant around the origin and this simplifies 
the problem drastically. In spite of the action complexity due to the realistic physical sit- 
uation involved, it will be always an easy way to find the quadratic expansion around the 
saddle color points and the resultant integration over the color variables is evaluated using 
the standard Gaussian quadrature. 

The critical point for the Hagedorn phase transition is determined in the midway of the 
interpolation between the two different analytical solutions for the small- and the large- A/A^, 
respectively. The Gaussian-like saddle points method is compared with the spectral density 
of color eigenvalues method (i.e. spectral density method, see for example Refs. 
and the exact numerical solution as well. It is found that the exact numerical solution 
fits precisely the results of the Gaussian-like saddle points method for the large thermal 
running parameter \/N% > Xq/N^ and furthermore it is found in a good agreement with the 
spectral density of color eigenvalues method. When the large- X/N^ is extrapolated to the 
small- A/iV^ < X /N^, the solution becomes slippery and when the X/N% reaches some critical 
value, it is deflected to increase. This solution is concave up and has a minimum at the critical 
point A ( - r/ . y m i Tr This point is the threshold point for the acceptable physical solution (II). 
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The minimum point for the extrapolation of the asymptotic large- A/Af? solution is actually 
the minimal threshold point. The extreme left hand side interval for the asymptotic large- 
\/N% solution is presumed to start from the threshold point A ( . // ^ m j n . The extrapolation 
down below this threshold point is unphysical and is declined as a solution. This means 
that beyond the threshold point, the asymptotic large- A/Af? solution (II) is deflected and 
changes its analyticity in order to match and satisfy the asymptotic small- A/A 7 ? solution (I). 
The solution (I) is found in agreement with the exact numerical results for the small A/A 7 ?. 
The asymptotic small- A/A^ solution is found be adopting standard approximations. The 
action is expanded to A-power expansion and then is evaluated using the group orthogonality 
over the full color range. The ultimate limit of the approximation validity is given by the 
point Ammax an d beyond this point the solution will be broken. The action structure is 
simplified drastically in the range A < A(/)max- For example, the quantum statistics can be 
approximated always to Maxwell-Boltzmann statistics and so on. However, the asymptotic 
small- and large-A/A 7 ? solutions (i.e. solutions (I) and (II), respectively) may split by a small 
additional constant. This constant is simply the approximation redundant in particular when 
both solutions are extrapolated far away from their asymptotic limits. The critical point 
Ao/A^ for the phase transition from the asymptotic small- A/A 7 !? to large-A/A 7 ? solutions is 
located near or above A > ^^^ m [ n but below A < A(/)max- It is roughly the midway 
interpolation between the small and large solutions. The advantage of the Gaussian-like 
saddle points method is the simplicity in the sense it can deal with a complicated physical 
problem such as the deconfinement phase transition in QCD. 

The QCD-phase transition is studied in the context of the color-singlet state of quark and 
gluon bags. It is found that the density of states for low-lying masses is a discrete spectrum. 
The low-lying density of states is determined by the known hadronic mass spectrum particles 
that are found experimentally and are available in the data book [45|]. The Hagedorn states 
appear just above the low-lying known mass spectrum. The mass spectrum for the Hagedorn 
states is continuous. The gas dominated by the low-lying mass spectrum particles undergoes 
a third order phase transition to a hadronic phase dominated by the Hagedorn states. When 
the system is thermally excited beyond the Hagedorn phase, the hadronic phase undergoes 
another phase transition to an explosive quark-gluon plasma. 

The scenario for the phase transition is depicted in Fig. ffTT]) . In the diluted nuclear 
matter when the temperature increases and reaches the critical one, the gas of the discrete 
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low-lying hadronic mass spectrum particles undergoes a third order phase transition to a gas 
of continuous Hagedorn states (i.e. the high-lying hadronic states). Furthermore, when the 
system is thermally excited above this temperature, the gas of Hagedorn states undergoes 
another first order or higher order phase transition to an explosive quark-gluon plasma. 
Moreover, the scenario for the phase transition can be extended to consider into account the 
multiple intermediate processes. An alternate scenario is depicted in Fig. (fl2|) . The gas of 
Hagedorn states undergoes a higher order phase transition to a gas of neutral colored bags 
or even non-singlet bags. These bags are not color-singlet states but carry neutral color 
charge and this charge is fixed by the effective color chemical potentials. Subsequently, the 
gas of neutral colored bags undergoes a higher order phase transition to a gas of colored 
bags. The non-singlet bags are not colored bags and not true deconfined colors but rather 
bags with specific internal color structures. In the colored bags, the quarks and gluons are 
still bounded in finite size blobs and they carry total color charge. The phase transition to 
the colored bags with conserved color charges might be associated with breaking the group 
symmetry SU(N C ) to {/(l)^ 0-1 . These colored bags become unstable and when the system 
is perturbed slightly and thermally excited, the system eventually undergoes another higher 
order phase transition to a real deconfined quark-gluon plasma. 

The possible consistent color- deconfinement and chiral restoration phase transition sce- 
narios are depicted in Fig. ffTBl . It is illustrated that the discrete low-lying hadronic mass 
spectrum undergoes a third order phase transition to the Hagedorn phase. The gas of contin- 
uous Hagedorn states with internal structures a > 2 undergoes a first order phase transition 
to an explosive quark-gluon plasma. However, it is possible in some alternate scenarios that 
the Hagedorn phase undergoes a higher order phase transition to the gas of non-singlet bags 
(i.e. with a < 2) and the subsequent gas of exotic states undergoes a higher order phase 
transition to an explosive quark-gluon plasma. The chiral phase transition for the gas of 
Hagedorn states with the internal structure a = 3 is found of the second order and it takes 
place utmost simultaneously the deconfinement first order phase transition to an explosive 
quark-gluon plasma. However, in the case that the chiral phase transition persists to take 
place in the gas of discrete low-lying hadronic mass spectrum states prior to the Hagedorn 
phase then the chiral phase transition will be a smooth cross-over one. In this case, it will be 
hard to reach the Hagedorn threshold production and subsequently the Hagedorn states will 
not be produced and no explosive quark-gluon plasma will be generated. The quark-gluon 
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plasma in this scenario can be found by the Gibbs construction. It is possible in other sce- 
narios that the the chiral phase transition coincides the Gross- Witten point and it will be a 
cross-over transition. It seems that the chiral restoration phase transition likely takes place 
utmost in the Hagedorn phase just below the deconfinement phase transition or at least 
coincides the Gross- Witten point. In some scenarios, it is likely that the chiral phase tran- 
sition coincides the Gross- Witten point while in other scenarios the chiral restoration phase 
transition is likely takes place simultaneously with the deconfinement phase transition to an 
explosive quark-gluon plasma. Indeed, the Hagedorn bag's internal structure a is found es- 
sential in the phase transition diagram. Furthermore, it is also possible the Hagedorn phase 
undergoes a higher order phase transition through multi-processes internal-structure phase 
transition while the chiral restoration phase transition is a smooth cross-over transition and 
takes place in the continuous Hagedorn phase. The order of the chiral phase transition 
becomes sharper and of a lower order as its point approaches the deconfinement point to an 
explosive quark-gluon plasma. The QCD phase transition diagram is proved very rich and 
non-trivial. 
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FIG. 1: (Color online) The saddle points solution for particles in the fundamental representation. 
The saddle points are distributed uniformly over the range \0i\ < ir for a parameter runs over the 
range A < Ao- The numerical solution shows that saddle points cease to exist when the thermal 
running parameter A exceeds the value = | in the large iV c limit. 
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FIG. 2: (Color online) The density of states for the color-singlet fireballs (i.e. Hagedorn states) 
versus the thermal running parameter with various color numbers N c . The high- lying solution 

^ c 

is calculated for a bag consisting of fundamental particles with a runs over the range > |. 
The results for the high-lying solution which is calculated using the method of the spectral density 
of color eigenvalues introduced by Brezin et. el. [8| in the limit of large color number N c and 
the results for the high-lying solution calculated using the Gaussian-like saddle points method 
are displayed together with the exact numerical solution. The high-lying solutions for the both 
methods and their extrapolations to < (77?) ^ are compared with the exact numerical 

solution. 
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FIG. 3: (Color online) The density of states for the color-singlet fireballs (i.e. Hagedorn states) 
versus the thermal running parameter for various color numbers N c . The low- lying density 
of states is calculated for bags consisting of fundamental particles and with a thermal running 
parameter runs over the range < |. The results of the density of states derived by the 
orthogonal expansion given by Eqs. ([56|) and ([57|) are displayed together with the exact density of 
states evaluated numerically. 
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FIG. 4: (Color online) The match between the asymptotic low-lying energy for jfe runs over the 

^ c 

range < yj^z} j anci the asymptotic high-lying energy solution for runs over the range 

■j^z > (l^) ] f° r a color-singlet bag consisting of only fundamental particles. The asymptotic 

low-lying solution is calculated using the A-power expansion and the group orthogonality over the 
entire color range \9\ < ir. a) The asymptotic high- lying energy solution is computed using the 
Gaussian-like saddle points method, b) The asymptotic high-lying energy solution is computed 
using the spectral density of color eigenvalues method ^]. 
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FIG. 5: (Color online) The high-lying density of states for the color-singlet bag consisting of only 
adjoint particles versus the thermal running parameter for various color numbers N c . The 
asymptotic high-lying energy solutions over the range -^j > (f ^p^ n = \ and their extrapolation 
to the small -hz are displayed and compared with the exact numerical one. The asymptotic high- 
lying energy solutions are computed using the Gaussian-like saddle points method and the spectral 
density of color eigenvalues. 
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FIG. 6: (Color online) The asymptotic low-lying energy solution for a bag of adjoint particles in the 
color singlet state in the limit N c — > oo versus the thermal running parameter -^j . The asymptotic 
low-lying energy solution versus is compared with the exact numerical one for various color 
numbers JV C =3, 4 and 5. The exact numerical solution is found to fit the asymptotic low-lying 
energy solution over the range < (f ^P^ n = \ < cr j^ ca ^ ■ The low- lying energy solution 
is found vanished as Mxa.^^^ In f -^Low-lying (^/"^c )) — * 0- The asymptotic high-lying energy 
solution limjvj^oo In f -^High-lying (V-^c )) vanishes to match the low-lying solution at the point 
of phase transition. 



103 



0.7 



0.6 



0.5 



0.4 



^"cn 0.3 
bfj 
O 



0.2 



0.1 



0.0 



Numeric: Exact numerical evaluation 
Asympt (Low- lying): Analytical approximation f./ t 

Log(Z)=-log(l-A x )«A x + A 2 x /2+ ... 




-0.1 



0.0 



Numeric: N c =2 
Numeric: N c =3 
Numeric: N c =4 
Numeric: N c =5 
Log(Z)=-log(l-A K ) 
Log(Z)=A x +A 2 N /2+A 3 N /3 
Log(Z)=A x +A 2 N /2 



0.1 



0.2 



0.:! 



0,1 



0.5 



A N =A/N C 



FIG. 7: (Color online) The density of states for the color-singlet fireballs (i.e. the Hagedorn states) 
versus the thermal running parameter with various color numbers N c . It is calculated for 
color-singlet bags consisting of only adjoint particles. The exact numerical solution is compared 
with the asymptotic low-lying energy solution with various analytic approximations for the small 
thermal running parameter < ^^p" 1 = ^ < cr j^ ca ^ ■ The exact numerical solution is found 
to fit the asymptotic analytic solution Z(X) = 1/(1 — t^t). 
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FIG. 8: (Color online) The high-lying density of states and its extrapolation for the color-singlet 
bag consisting of fundamental and adjoint particles with the same thermal running parameter 
versus the thermal running parameter with various color numbers N c . The high-lying energy 
solution is the asymptotic solution for the thermal running parameter runs over the range > 
CT jV^ ca ^ — (P re ^ r ^ ca ^ > (/ ~^ mm , The asymptotic high- lying energy solutions are computed 
using the Gaussian-like saddle points method and the spectral density of color eigen-values and 
they are compared with the exact numerical one. 
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FIG. 9: (Color online) The asymptotic low-lying energy solution for color-singlet bags consisting 
of fundamental and adjoint particles with the same thermal running parameter for the both 
fundamental and adjoint particles versus with various color numbers N c . The exact numerical 
solution is compared with the analytic approximation of the low-lying energy solution for the 
thermal running parameter runs over the range < ^. It is found that the exact numerical 
solution fits the approximate analytical one for small values of When the thermal running 
parameter exceeds the (pre-)critical one > (j^)fp re ) cr itical' ^ e exac * numer i ca l solution 
starts to deviate slightly from the approximate analytical solution and this deviation becomes 
noticeable beyond the point A critjcal _ ^_A^ = 0.34. 
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FIG. 10: (Color online) The asymptotic low-lying and high-lying energy solutions and their ex- 
trapolations for the color-singlet bag of fundamental and adjoint particles with the same thermal 
running parameter versus for various color numbers N c . The exact numerical solution fits the 
asymptotic low-lying energy solution over the range < cr j^ ca ^ and fits the high-lying energy 
solution over the range -j^ > critical . The asymptotic low- lying and high- lying energy solutions 
split by a small redundant constant at the critical point of the phase transition. This redundant 
small constant emerges due to the kind of approximations are considered. 
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FIG. 11: The phase transition scenario in the context of the present model. At low temperature, the 
diluted nuclear matter is dominated by the low-lying hadronic mass spectrum such pions, nucleons 
etc ■ ■ ■. At the Hagedorn critical temperature, the system undergoes a third order phase transition 
to system dominated by the highly excited hadronic states known as Hagedorn states or fireballs. 
When the system is thermally excited to higher temperature, the hadronic phase undergoes a first 
order phase transition to an explosive quark-gluon plasma. 
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FIG. 12: Same as FigQT] but with an alternate scenario for the phase transition. The hadronic 
phase dominated by Hagedorn states undergoes a higher order phase transition to another phase 
dominated by highly excited neutral color bound quark-gluon bags (i.e. bound state but with a 
color non-singlet state). At higher temperature, the gas of colored quark and gluon bags undergoes 
a higher order phase transition to an explosive real deconfined quark-gluon plasma. In this scenario, 
the phase transition to an explosive quark-gluon plasma takes place only through more complicated 
multi-processes. 
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FIG. 13: The consistent chiral and deconfinement phase transition scenario. 



110 



